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Abstract 

In this report we review recent developments in perturbation theory methods for gauge 
theories. We present techniques and results that are useful in the calculation of cross sections 
for processes with many final state partons which have applications in the study of multi-jet 
phenomena in high-energy Colliders. 
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1 Introduction 



In high energy collisions among hadrons and/or leptons the production of final states with a large 
number of energetic, widely separated partons gives rise to events with many jets in the final state 2 . 
In many cases these multi-jet events offer a potentially important probe on new physics |29| . e.g. in 
the case of the sequential decays of new heavy particles, such as a Higgs decaying to four jets through 
real W/Z pairs, or such as a pair of heavy gluinos decaying into a multi-jet system through a chain- 
decay of the various unstable supersymmetric particles. The possibility of using these observables 
to identify new phenomena relies on our capability to predict the production rates and features 
of the standard multi-jet production mechanisms which often provide a significant background to 
these discovery channels. 

Monte-Carlo techniques exist to describe processes with many partons in the final state through 
a branching process driven by the leading logarithmic approximation to the multiple emission 
probabilities j^H]- This approach provides a scheme in which the number of final state particles 
is not fixed, and can be as large as allowed by the relative branching probabilities. Most of the 
emitted particles will be either soft or collinear to the leading ones involved in a primary scattering 
process (say gg — > gg), because these configurations are enhanced by the dynamics. Inclusive 
energy measurements, such as calorimetric detection of jets, are insensitive in first order to the 
precise details of the structure of a collinear shower and, in the case just mentioned of a gg —>■ gg 
scattering, will usually only detect two jets at large angle. Events with multi-jets are generated in 
the Monte-Carlo approach whenever a branching with large relative transverse momentum takes 
place. In this case a new branch of the partonic shower will arise and will independently evolve as 
a secondary jet. However, while the branching probabilities properly describe the parton evolution 
within a jet in the leading log approximation, this approximation does not properly describe the 
emission of partons at large relative transverse momentum. For these processes, therefore, a full 
calculation of the matrix elements for the hard process involving the many leading partons is 
required. 

The calculation of these processes is made particularly difficult by the large number of Feynman 
diagrams which appear in the perturbative expansion. As an example in Tabled we have collected 
from Ref. [53] the number of diagrams contributing to the process gluon-gluon — > n-gluons. 

The structure of the non-abelian vertices, furthermore, leads to an almost uncontrollable inflation 
in the number of terms which are generated, and very soon standard techniques of numerical 
evaluation or algebraic symbolic manipulation become useless. Significant simplifications in these 
calculations have been achieved in recent years thanks to the use of new simple representations 
for vector polarizations, a better organization of the diagrammatic expansion which fully exploits 
the properties of gauge invariance, the discovery of recursive relations which connect amplitudes 
with n + 1 partons to amplitudes with n ones and the use of supersymmetric Ward identities 
to relate gluonic and quark-gluon amplitudes. In spite of these advances, the results of these 
calculations are still often very complicated and sometimes of limited use, even numerically, for 

2 For experimental analyses of multijet production in e + e~ collisions sec, for example, [201 IS1 IHSI - For production 
of multijets in pp collisions, see jHHl El 1231 IH8I I91j . For associated production of weak gauge bosons and jets in pp 
collisions, see IH3 EDI 1221- 
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systematic analysis of their phenomenological implications. In addition to the development of these 
tools for the calculation of exact matrix elements, effort has therefore also been put into finding 
proper approximations which reliably simulate the exact solutions in the relevant regions of the 
multi-particle phase-space and which are sufficiently simple to be handled analytically and fast to 
evaluate numerically. 



n 


2 


3 


4 


5 


6 


7 


8 


# of diagrams 


4 


25 


220 


2485 


34300 


559405 


10525900 



Table 1: The number of Feynman diagrams contributing to the scattering process gg — > n g . 

In this Report we collect and review these recent developments for the calculation of multi-parton 
matrix elements in non-abelian gauge theories. For examples of how these matrix elements can be 
used to obtain cross sections for processes in high energy colliders see EHLQ [29J and references 
contained within. 

In Section 2 we describe the helicity-amplitude technique and introduce explicit parametrizations 
of the polarization vectors in terms of massless spinors. To reach a wide an audience as possible we 
have chosen not to use the Weyl - van der Waarden formalism preferred by some researchers, see 
for example Ref.|10j. 

In Section 3 we introduce an alternative to the standard Feynman diagram expansion, based 
on the equivalence between the massless sector of a string theory and a Yang-Mills theory. This 
expansion groups together subsets of Feynman diagrams for a given process in a gauge invariant 
way. These subsets are easier to evaluate than the complete set and different gauges can be used for 
each subset so as to maximize the simplifications induced by a proper choice of gauge. Furthermore, 
different subsets of diagrams are related to one another through symmetry properties or algebraic 
relations and can be obtained without further effort from the knowledge of a small number of building 
blocks. This expansion can be extended to arbitrary processes involving particles in representations 
other than the adjoint, and in this Section we construct this generalization. 

Section 4 describes the use of Supersymmetry Ward identities to relate amplitudes with parti- 
cles of different statistics. These relations are useful even when dealing with non-supersymmetric 
theories because in many cases the additional supersymmetric degrees of freedom decouple from 
the processes of interest. In addition, if the energy of the scattering process is large with respect to 
the mass splittings within supersymmetry multiplets, these relations can be used to easily calculate 
the matrix elements for the production of supersymmetric particles. 

In Section 5 we illustrate the use of these tools with the explicit calculation of matrix elements 
for processes with four and five partons, and give results for the scattering of six gluons and four 
gluons plus a quark-antiquark pair. We hope this Section is useful for the reader who wants to 
familiarize himself with the details of how these calculations are performed. 

In Section 6 we prove various factorization properties using a string-theoretic approach, which 
provides a compact way to represent multi-parton amplitudes. The results contained in this Section 
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are useful for a better understanding of the structure of multi-parton amplitudes in gauge theories. 

Section 7 introduces the Berends-Giele recursion relations, which allow to calculate matrix el- 
ements in a recursive fashion, providing an algebraic algorithm which can be efficiently used for 
numerical evaluation of higher order processes. 

In Section 8 we collect some explicit results concerning matrix elements for processes with an 
arbitrary number of particles. These expressions hold for amplitudes with a simple helicity structure, 
and whose properties are fully determined by their behaviour at the collinear and infrared poles. 
These results help understanding and extending known coherence properties of the soft radiation 
in non-abelian gauge theories, as will be discussed. 

In Section 9 we show how to use these techniques in the case of gauge groups which are the 
product of different groups, and how to calculate in presence of massive gauge bosons from a 
spontaneously broken gauge theory. As an application, we collect the known matrix elements for 
the processes involving a massive gauge boson produced in association with two gluons from the 
scattering of a quark-antiquark pair. 

Section 10 describes the approximation techniques mentioned above. We review different ap- 
proaches that have been proposed and illustrate their use for processes involving n-gluons or n-gluons 
plus a quark-antiquark pair. 

Finally, we collect in five Appendices various definitions, conventions and results which are useful 
in performing explicitly analytic calculations or numerical evaluations of matrix elements. 
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2 Helicity Amplitudes 



The use of helicity amplitudes for the calculation of multi-parton scattering in the high-energy 
(massless) limit was pioneered in papers by J.D. Bjorken and M. Chen [21], and by 0. Reading- 
Henry |H3j, and later further developed and fully exploited by the Calkul Collaboration in a classical 
set of papers (221 CHI- The application of this technique to QED processes is extensively reviewed 
in the book by Gastmans and Wu • 

According to this approach, one calculates matrix elements with external states having a given 
assigned helicity Since different helicity configurations do not interfere, to obtain the full cross 
section it is sufficient to sum incoherently the squares of all of the possible helicity amplitudes 
which can contribute to the process. The advantage over more standard techniques is that by 
choosing a definite helicity configuration one can exploit gauge invariance and select an explicit 
representation for the polarization vectors which will simplify the calculation. 

Since the polarization vectors always enter in an amplitude contracted with a gamma matrix in 
QED processes, the Calkul group found it useful to introduce a representation in terms of the two 
momenta (q, q') of one of the pairs of external charged fermions in the process : 

e±(p) = e±(p) • Y = N[pq'q(l ± 7s ) - q'qp(l T 7s) T 2(g • </)P7 5 ], (2.1) 

where N is a normalization factor, 

N = [lQ(q-q')(p-q)(p-q')]-^ (2-2) 

With this choice for e , many of the terms which appear in the diagrammatic expansion simply 
vanish. Because of gauge invariance, all of the terms generated by the third piece in Eq. ()2.1|) . 
proportional to p, will sum up to zero. Furthermore, helicity conservation along a fermionic line 
guarantees that at least one of the two remaining terms, containing orthogonal chiral projections, 
will vanish. Finally, if the photon happens to be attached to an external fermion whose momentum 
is one of the reference momenta used to define the photon polarization, then this diagram will also 
vanish provided the helicities match (see the e + e~ annihilation case in the Appendix for an explicit 
example). 

If the set of diagrams contributing to the given matrix element can be split into the sum of gauge 
invariant subsets, we can choose different reference momenta q, q 1 for different subsets, provided we 
keep track of the relative phase which can appear in the polarizations when they are referred to 
different g's. As an example of a process in which this splitting is possible, we indicate e + e _ — > 
/i + yU~7, and refer to the previously quoted papers for the explicit calculation. 

While this technique turns out to be extremely useful for pure QED calculations, the complexity 
of a non-abelian theory calls for something even simpler. In the non-abelian theory, in fact, the 
proliferation of diagrams is such that the bookkeeping of the different phases becomes very complex, 
and the existence of processes without external fermions calls for a different choice of reference 
momenta to achieve the desired simplification. An improved version of the Calkul representation, 
which is more apt to use in non-abelian theories, was introduced by Xu, Zhang and Chang in Ref. [§2] 
and, independently, in Ref. I n this improved version, a vector polarization is expressed in 
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terms of massless spinors and just one reference momentum. Here in the following we will give 
a simple derivation of this result making use of Supersymmetry [7j. We will always work in four 
space-time dimensions, but the construction could be extended in principle to higher dimensions, 
and possibly to non-integer dimensions as well. 

To start with, we will set our notation and will present some definitions concerning the spinor 
algebra that will be extensively used in the following. For additional details and properties, see the 
Appendix. 

Let ip(p) be a massless four-dimensional Dirac spinor, i.e. : 

p i/}(p) = p ■ Jip(p) = p 2 = 0. (2.3) 
We define the two helicity states of i/j(p) by the two chiral projections: 

MP) = £(1±75MP) = MpT, (2-4) 

the last identity being just a conventional choice of relative phase between opposite helicity spinors 
fixed by the properties under charge conjugation ( c ): 

^(p) c = &/>(?)• , C^C- 1 = v (2.5) 

Following [52], we introduce the following notation: 

\p±) = V±(p) (p±| = MP) (2.6) 

(pq) = (p-\q+) = ip_(p)ip + (q) [pq] = (p + \q-) = ^+(p)^_(g). (2.7) 

The spinors are normalized as follows: 

(p±| 7M |p±) = 2jy (2.8) 
From the properties of the Dirac algebra, it is straightforward to prove the following useful identities: 

(p+\q+) = (p-\q~) = {PP} = [PP] = (2.9) 



(pq) = -(qp), [pq] = -[qp] (2.10) 
\(pq)\ 2 = 2(p-q) (2.11) 

(A + 1 7m \B+) (C-\Y \D-) = 2[AD](CB). (2.12) 

We now turn to the description of four- dimensional massless vectors. In four dimensions the 
physical Hilbert space of a massless vector is isomorphic to the physical Hilbert space of a massless 
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spinor (up to a Z 2 transformation), since they both lie in one-dimensional representations of 50(2), 
the little group of S0(3, 1). This isomorphism is realized through a linear transformation which 
relates like-helicity vectors and fermions: 



e 



(p) = Au + (p) lfl v, (2-13) 



e~{p) = (e+(p))' (2.14) 

where ez(p) is the polarization vector of an outgoing (i.e. positive-energy-) massless vector of 
momentum p, u + (p) is a massless spinor as defined above, v a is an a priori arbitrary Dirac spinor 
and A is a normalization constant, needed to enforce the usual normalization conditions: 

e+(p) • e+(p) = 0, e + (p) - e~(p) = -1. (2.15) 

In this isomorphism, the gauge invariance associated with the massless vector can be parametrized 
by the arbitrariness in the choice of the spinor v. Although this parametrization does not exhaust 
all the possible gauge choices, nevertheless it will turn out to be particularly useful in the following. 
It is easy to check that by properly choosing the gauge we can always select a spinor v(k) to be 
used in (|2.13|) that satisfies the following properties: 

kv(k) = k ■ iv(k) = 0, (2.16) 

k 2 = 0, k-p^O. (2.17) 

We will refer to the arbitrary k as to the reference momentum. Therefore we can always write, for 
a proper gauge choice: 

e+(p,k) = A(p+\^\k+) (2.18) 

e~(p,k) = A* {p-^k-) (2.19) 

The normalization A has to be chosen to give unit norm to the polarization. Using Eq. (j2.12|) 
we easily obtain: 

e+(p,k) -e-(p,k) = -2\A\ 2 p ■ k (2.20) 



and thus: 



e+(p,k) = e W> ^ + (2.21) 



V2(kp) 



(2.22) 
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where <fi(p, k) is a phase which a priori depends on the vector momentum p, and on the reference 
momentum k. If we set this phase to zero, it is easy to show that that the change in the polarization 
vector caused by a change in the reference momentum is given by: 

^(p,k)^4(p^)-y^ 1 j^ j p ll . (2.23) 

Note that for this choice of <f>(p, k) the a priori phase factor in front of e+(p, k') is equal to unity 
in this equation. A similar result holds for the negative helicity vectors. Therefore the choice of 
polarization vectors used through out this review is 

4( P ,k) = ± fr*' ^ (2.24) 

Using this representation, ([2.24)1 . for the polarization vectors in the calculation of a given ampli- 
tude, we can choose not only a different reference momentum k for each polarization vector in the 
process, but we can also choose different reference momenta for each gauge invariant part of the full 
amplitude, without having to worry about relative phases. This property will be used extensively in 
the following applications, where we will decompose each amplitude into a sum over gauge invariant 
components. 

A proper assignment of reference momenta to the different external vectors will result in signif- 
icant simplifications. As an example, by using Eqs. ([2. 12)1 . ([2.9)1 one can easily prove the following 
identities: 

e+(p, k) ■ e+(p', k) = e + (p, k) ■ e~{k, k') = (2.25) 

These identities suggest that it is convenient to choose the reference momenta of like-helicity vectors 
to be the same and to coincide with the external momenta of some of the vectors with the opposite 
helicity. 

The representation (|2.24|) for the polarizations is also particularly helpful when calculating 
processes with external fermions in addition to the vectors. The polarization vectors contract 
with the gamma matrices in the following way: 

6 ± (p,fc)- 7 = ± — ^ (\ PT )(k T \ + \k±)(p±\). (2.26) 
{kT\P±} 

An explicit example of the use of these formulas for the simple case of e + e~ annihilation into two 
photons is given in the Appendix. 

As a final comment, we add that the gauges generated by this choice of polarization vectors are 
equivalent to axial gauges. In fact it is straightforward to prove on the basis of the identities given 
here and in the Appendix, that: 

( e tT = e+(p,k)e-(p,k) + e-(p,k)e+(p,k) = -g, u + P^ +P"^ . ( 2 .27) 

pol " 

Because of this reason, we will expect these gauges to make calculations particularly simple when 
studying matrix elements in the eikonal approximation. 

The representation of polarization vectors in terms of spinors has been generalized to the case 
of massive particles of spin 1/2, 1 and 3/2 in Ref . [82) . 
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3 The Color Form Factors 



3.1 Gluonic Amplitudes: Duality and Gauge Invariance 

In perturbative QCD the calculation of multi-gluon scattering amplitudes, even at tree level, is 
very challenging. The number of diagrams describing a given process grows very quickly, and 
the redundancy due to the gauge invariance leads to a rapid proliferation of terms. One way to 
simplify these calculations is to divide all of the diagrams contributing to a given matrix element into 
subsets of diagrams which are independently gauge invariant under redefinition of the polarizations: 
e i (Pi) ~> e ii.Pi) + a iPi, with the CKj's being arbitrary functions. It might then be possible to 
choose different gauges for these different subsets in such a way as to simplify the calculation as 
much as possible. By using the polarization vectors introduced in the previous Section, different 
gauge choices will not change the relative phases between the different gauge invariant pieces, thus 
contributing to a further simplification. 

The issue then is to find a systematic way of dividing processes into gauge invariant components. 
In this Section we will provide such a criterion based on the work initiated in jOHl ITUl loT?] . This 
criterion can be applied to any gauge theory: here for simplicity we will refer to simple unitary 
groups SU(N), but the techniques introduced can be easily extended to more general cases, such 
as products of groups, as will be shown in a later Section. 

A very complete study of the relation between gauge invariance and color structures in the 
context of the large- N limit jlH] of QCD and the loop expansion was presented by Cvitanovic and 
collaborators in Ref.j2Zj. Some of the results presented here do overlap with theirs. 

For the sake of reference, we will often refer to the Yang-Mills gauge bosons as to gluons. As 
we will prove in what follows, it turns out to be useful to consider the space of color configurations 
for the given scattering process. If we expand the amplitude with respect to an orthogonal basis in 
this space, this expansion is guaranteed to be gauge invariant. Therefore there are many different 
ways of breaking up the amplitude into gauge invariant components. A particular choice which can 
be singled out for its prompt physical interpretation and for its many important properties is to 
insure that these gauge invariant components be invariant under cyclic permutations of the external 
gluons. Consider an SU(N) Yang-Mills theory; then at tree level in perturbation theory any vector 
particle scattering amplitude, with colors 0,1, 0,2 ■ ■ ■ O n , external momenta pi,p 2 - ■ - Pn and helicities 
ei, e 2 . . . e n , can be written as 

M n = Yl tr{X ai X a2 ...X an )m{p 1 ,ei;p 2 ,e 2 ;---;p n ,e n ), (3.1) 

{1,2,. ..,n}> 

where the sum with the prime, D{i,2,...,n}'> is over an i n ~ 1)' non-cyclic permutations of 1, 2, . . . , n 
and the A's are the matrices of the symmetry group in the fundamental representation, which we 
choose to normalize as follows 3 : 

[A a , \ b ] = if abc \ c , tr(\ a \ b ) = S ab . (3.2) 

3 This normalization of the A matrices differs from the usual one by a y/2, which we explicitly add to the Feynman 
rules (see Appendix C): this choice is purely conventional, and just simplifies the bookkeeping of factors of 2 in the 
calculations. 
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The color structures given by the traces of A matrices do not provide a complete basis for the 
possible color configurations of n gluons, but nevertheless they are sufficient to describe the tree- 
level scattering of n-gluons, as we will show below. It should also be pointed out that the color 
structures used in equation (|3.1|) are only orthogonal at the leading order in the expansion in 
powers of N; if {a} and {b} are two permutations of the gluon color indices we have in fact (see 
the Appendix): 

tr(\ ai \ a2 . . . A a ") [tr(\ bl \ b2 . . . X b ")]* = N n ~ 2 (N 2 - 1) ( 5 {a}{b} + 0{N- 2 ) ), (3.3) 

ai =l,jV 2 -l 

where the £{a}{6} is equal to 1 if and only if the two permutations are the same (up to cyclic re- 
orderings); this partial orthogonality, nevertheless, is clearly still sufficient to guarantee the gauge 
invariance of the expansion, which must hold order by order in 1/N. For a different choice of base in 
the color space, which is exactly orthogonal, see the alternative approach developed by Zeppenfeld 
in Ref.l 



The proof of Equation (|3.1|) is very simple if one uses the relations ()3.2j) : in any tree level Feyn- 
man diagram, replace the color structure function at some vertex using f abc = —i tr(X a X b X c — X c X b X a ) 
Now each leg attached to this vertex has a A matrix associated with it. At the other end of each of 
these legs there is either another vertex or this is an external leg. If there is another vertex, use the 
A associated with this internal leg to write the color structure of this vertex f C( i e X c as — % [X d , X e ] . 
Continue this processes until all vertices have been treated in this manner. Then this Feynman 
diagram has been placed in the form of Equation (|3.1|) . Repeating this procedure for all Feynman 
diagrams for a given process completes the proof. 

The sub-amplitudes m(l, 2, . . . , n) = m(pi, ei,p2, ti\ ■ ■ ■ Pn, Cn) of Equation (j3.1|) are by con- 
struction independent of the color indices and satisfy a number of important properties and rela- 
tionships: 

1. m(l, 2, . . . , n) is gauge invariant. 

2. m(l, 2, . . . , n) is invariant under cyclic permutations of 1, 2, . . . , n 

3. m(n,n — 1, . . . , 1) = (— l) n m(l,2, . . . ,n) 

4. The Dual Ward Identity: 

m(l, 2, 3, . . . , n) + m(2, 1, 3, . . . , n) + m(2, 3, 1, . . . , n) 

+ ■•■ + m(2,3,...,l,n) = 0. (3.4) 

5. Factorization of m(l, 2, • ■ • , n) on multi-gluon poles. 
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6. Incoherence to leading order in the number of colors: 

£ \M n \ 2 = N n ~\N 2 - 1) £ {|m(l,2,---,n)| 2 + 0(iV- 2 )}. (3.5) 

colors {l,2,...,n}' 

This set of properties for the sub- amplitudes we will refer to as duality and the expansion in 
terms of these dual sub-amplitudes the dual expansion. Properties (1) and (2) can be seen directly 
from the properties of linear independence (to the leading order in N, and for arbitrary N) and 
invariance under cyclic permutations of tr (A 1 A 2 . . . A n ). Whereas (3) and (4) follow by studying the 
sum of Feynman diagrams which contribute to each sub- amplitude. The sum of Feynman diagrams 
which enter into the Dual Ward Identity is such that each diagram is paired with another with 
opposite sign so that the combination contained in Equation ()3.4|) trivially vanishes. Property (5) 
will be discussed in detail in Section 6 and the incoherence to leading order in the number of colors 
(6) was obtained above and follows from the color algebra of the SU(N) gauge group. 

To the string theorist this expansion and the duality properties (1) to (6), see jJU], are quite 
familar since the string amplitude, in the zero slope limit, reproduces the Yang-Mills amplitude on 
mass shell |85j . Each sub-amplitude then corresponds to the zero slope limit of a string diagram, and 
the sub-amplitude can be obtained by using the usual Koba-Nielsen formula [S3]. Kawai, Lewellen 
and Tye, Ref.[22], have derived a relationship between the closed string tree amplitudes and the 
open string tree amplitudes which allows this connection to be explicitly extended to the heterotic 
string as well as to the closed bosonic and the type II superstring. The traces of A matrices are just 
the Chan-Paton factors |83j . For the string amplitude the properties (1) through (6) are satisfied 
even before the zero slope limit is taken, and in particular Equation ()3.4|) holds as a Ward identity 
for correlation functions of products of two-dimensional conformal fields. We will see later on in 
this Section and also in Section 6 on factorization properties, useful examples of how to use the 
string representation to derive various properties of the Yang-Mills amplitudes. 

Which diagrams contribute to a given sub-amplitude and with which coefficients they enter can 
be determined by the procedure developed earlier in this Section for re-writing the color factors. 
It is however helpful to think in terms of string diagrams, and to realize that the contributing 
Feynman diagrams can just be obtained by pinching in all possible ways on multi-particle poles the 
string diagram itself (see for example Figure 1). 

The relationship with the string diagrams, the possibility of choosing an ad hoc gauge and the 
simple factorization properties that the dual sub-amplitudes must satisfy, suggest that a Yang- Mills 
amplitude expressed as in Equation (J3.1|) will assume a particularly simple form. That this is in 
fact the case will be shown in Section 5 , where we will consider some explicit examples. 

The gauge invariance and properties under cyclic and reverse permutations (12 . . . n) — > {n . . . 21) 
allow the calculation of far fewer than the (n — 1)! sub-amplitudes that appear in the dual expansion. 
In fact the number of sub-amplitudes that are needed is just the number of different orderings 
of positive and negative helicities around a circle. Of course some of the sub-amplitudes vanish 
because of the partial helicity conservation of tree level Yang-Mills and others are simply related 
to one another through the properties (2) through (4). Kleiss and Kuijf in Ref. (22] have given a 
detailed, general accounting of the minimum number of independent gluonic subamplitudes that 
are needed for the n-gluon scattering. Their results is that (n — 2)! subamplitudes are independent. 
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3.2 Quark-Gluon Amplitudes 



In this Section we will extend the color representation introduced above to processes with fermions 
in the fundamental representation of the gauge group SU (N) [27J EH EH] • We will aim at a repre- 
sentation which satisfies the properties of gauge invariance and factorization to the leading order in 
l/iVc, Eq. (j3.5j) . As before, we will refer to the gauge bosons as to gluons, and to the fermions as to 
quarks. We will start from processes with one quark-antiquark pair, and for the time being we will 
only consider tree level diagrams. 

The color structure of diagrams where all of the gluons are emitted directly from the quark line 4 
(see Figure 2) is obtained in a straightforward way directly from the Feynman rules: 

(A ai A a2 ...A a % (3.6) 

where are the color indices of the qq pair and ai v .. jn are the color indices of the gluons in 
the order they are emitted. In order to analyze diagrams with gluons coupling to each other, let 
us consider the case in which just one gluon is emitted from the quark and develops into a tree 
(Figure 3). We can factorize the color structure of the diagram into the color coefficient of the 
qqg vertex, namely (A A )jj, and the color structure of the remaining gluon tree. By using the dual 
representation, we can express this as the sum over traces of permutations of A matrices. As a 
result, we will obtain that the color coefficient of this diagram is given by a sum over permutations 
of the following expression: 

V(A A ) 4j tr[A A (A ai A a2 ...A a ™)] = (A ai A a2 . . . A a % - — <%tr(A ai A a2 . . . A a "). (3.7) 

A N 

This identity follows from the following property of the A matrices: 

n 2 -i y 

a=l JV 
with the normalization given in Equation (|3.2j) . 

The term proportional to A^ _1 corresponds to the subtraction of the trace of the U(N) group in 
which SU(N) is embedded. This trace couples to the quarks but commutes with SU(N) itself, and 
then it doesn't couple to the gluons. As such it must disappear after the sum over permutations. 
That this is in fact the case, can be easily checked. Since an arbitrary diagram can be factorized 5 
into diagrams of the QED type and diagrams with a tree evolution initiated by a single gluon, we 
conclude that any diagram with a qq pair can be decomposed in terms of the n\ permutations of the 
color structure given in Eq. (J3.6|) . Notice that all terms of the form 5ijtr(X ai A" 2 . . . A a ") do cancel 
(at tree level). 

By repeated use of the factorization properties of the color coefficients one easily arrives at the 
general representation in terms of which it is possible to decompose diagrams with more than one 

4 We will denote these diagrams as QED-type |27|. 

5 Again, the factorization we are referring to here is just a factorization of the color structure, and not of the full 
amplitude, since because of gauge invariance factorization cannot be applied on a diagram-by-diagram basis. 
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quark pair: 



A(K},{a}) = ^(A*...A^) 4l(n (A^+^..A^) faaa ...(A^+ 1 ...A°-) imata (3.9) 

Here m is the number of qq pairs present in the diagram, n is number of gluons emitted and 
the indices n\, . . . , n m _i (with 1 < n, < n) correspond to an arbitrary partition of an arbitrary 
permutation of the n gluon indices. A product of zero A matrices has to be interpreted as a Kronecker 
delta. The indices ii, . . . ,i m are the color indices of the quarks and the indices «i, . . . , a m are the 
color indices of the antiquarks. By convention all of the particles are outgoing, so each external 
quark is connected by a fermionic line to an external antiquark. When we want to indicate that a 
quark with color index and an antiquark with color index ol\. are in fact connected by a fermionic 
line, we identify the index a k with the index i k . Therefore the string {a} = (a±, . . . ,a m ) is a 
generic permutation of the string {%} = (ii, . . . ,t m ). The power p is determined by the number of 
correspondences between the string {a} and the string {t}, i.e. by the number of times a k = ik- 
If {«} = {2}, then p — m — 1. Contrarily to the process with only one qq pair, in which terms 
with p 7^ vanish after gauge invariant quantities are formed, here the terms with p 7^ do not 
vanish. The reason for this fact is that while one qq pair cannot couple via the £/(AQ-trace to a 
set of gluons, the U (AQ-trace can connect two qq pairs, and then it has to be explicitly subtracted 
if the gauge group we want is just SU(N). These subtraction terms are exactly given by the color 
structures in Eq. (j3.9|) proportional to A^~ p , (p > 0). The negative powers of A" are a consequence 
of the coupling between a quark and a £7(AQ-trace, which according to the normalization chosen in 
Ea.(l3~2D is given by l/\Hsf. 

To give an example, in the case of two quark pairs and two gluons the possible color structures 
are the following: 

(X a X b )i 1 i 2 5i 2 i 1 , (X a )i 1 T 2 (X b )i 2 t 1 , 5 il i 2 (X a X b )i 2 i 1J (3.10) 



^r(A A Si 2l2 , 7\r(A )iiii(A )i 2 i2 j at^«i«i (A X )i 2 i 2 , (3.11) 

where a and b represent the color indices of the two gluons, and the six additional color structures 
with a and b interchanged have been omitted. 

The representation given in Eq. (J3.9|) has the simple physical description which we will now 
illustrate. 

To start with, let us consider the color structure of an amplitude with quarks only, at tree 
level. As before, we will take all the particles as outgoing, and will assign indices (i 1 ,...,i m ) 
to the quarks and indices (ii, . . . ,i m ) to the antiquarks. It is understood that the quark i k is 
continuously connected through a fermionic line to the antiquark t^, for each 1 < k < m. Helicity 
will be conserved along this quark line, as well as flavor, since only gluons can be emitted. We 
can furthermore assume all quarks to be of different flavor, the case with identical quarks being 
similar but more confusing. It is easy to verify that the color functions accompanying each diagram 
contributing to this scattering process can be decomposed in terms of the following color structures: 

(-l) p 

D ({ a }) = <*iiaAa 2 • • • 5 i m a m (3-12) 
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where {a} = (a%, . . . , a m ) is a permutation of {?} = (Ti, . . . ,T m ). Each color structure D({a}) 
defines a color-flow pattern inside the diagram. As before, p is the number of [/(AQ-trace gluons 
that can appear in a given color-flow configuration, and whose contribution has to be subtracted. 
In this case, p is the number of 5i k ^ k appearing in the product, with the exception that when all 
of the delta functions connect quark pairs that belong to the same fermionic line (i.e. {a} = {i}) 
then p — m — 1. 

In Figure 4 the case m = 2 is shown, with the possible color factors given by (see Eq. (|3.8)0 : 

^UsAaJi > jy^tiJi'Wa (3.13) 

In Figure 4 the propagation of the {/(iV)-trace gluon is represented by the dashed line between the 
two quark lines. 

The color structure in Equation (J3.9|) has then a very simple physical interpretation. In fact it 
corresponds to the emission of the gluons off the color-flow lines defined by the functions D({a}). 
Each function D({a}) defines a net of color flows, as shown in Figure 5 for the case m = 2. Each of 
these color-flow lines, specified by a pair of indices (ik,Jk>), ac ts as a sort of antenna, that radiates 
gluons with an associated color factor (X...X) ik j, (see Figure 5). This color factor is the one 
appearing in the QED-type diagrams, i. e. diagrams in which all the gauge bosons are emitted from 
the fermionic line and no three- or four- vector vertices are present. Equation ()3.9j) shows that even 
graphs with non-abelian vertices can be decomposed as sums of QED-like diagrams. 

Given a helicity configuration for the external states the matrix element for m quark-pair plus 
n gluon scattering can then be expressed as: 

M m , n = J2 A p(MAa})rh p {niUa} (q,h). (3.14) 

where the sum is over the permutations of (ji, . . . , j m ). A p ({nj}, {a}) are the color factors appearing 
in Eq. ()3.9|) : they depend upon the partition and permutation {n«} of the gluon indices and upon 
the antenna pattern determined by the permutation of indices {a}. We introduced the subscript 
p to remind that depending on the permutation {a} the color factor will be proportional to a 
given power N~ p . The sub- amplitudes m p (q, h) multiplying a given color factor are functions of 
the momenta q and helicities h of the external particles. These sub-amplitudes are obtained by 
summing contributions from various different Feynman diagrams. If some of the external states are 
in a given color configuration, for example in a color singlet, the amplitude can be easily obtained 
by contracting Eq. (j3.14j) with the proper projector. 

To the leading order in N only the terms in the sum with ^ jj. will contribute, and the sum 
over colors of the amplitude squared will be the sum of the squares of the functions wt.?^ | a |, the 
interferences being suppressed by negative powers of N, as can be easily checked using Eq. ()3.8|) : 

El-^i 2 = iv m+ " E Kkw(^)I 2 (3-15) 

col {rii}{a} 

The 'hat' restricts the sum to the permutations {a} with a fe ^ j k for all fc's. 
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Each sub-amplitude r a \(q, h) is invariant under gauge transformations of the gluon polar- 
izations — > + /3p* . To prove this it is sufficient the orthogonality, to the leading order in N, of 
the color factors. We will now prove this fact. 

Let Sm^, r a \(q, h) be the gauge variation of a given sub-amplitude. Let {n,}, {a} be a given 
partition and a given permutation of quark and gluon indices chosen in such a way that p = 0. 
Then the following identity follows: 

=E A {?},w(ff^)^« = N n+m 5m {niUet} ( q ,h)+O(l/N 2 ). (3.16) 

col 

This shows that all the sub- amplitudes labelled by p = are gauge- invariant, since gauge- invariance 
does not depend on N and variations of 0(1/ N 2 ) cannot cancel the leading piece. We can now 
select all of the sub-amplitudes corresponding to p = 1, and repeat the same construction to show 
that they are gauge invariant too. In this way one can continue until p = m — 1 is reached, thus 
proving that each sub- amplitude is in fact gauge invariant. 

This gauge invariance is particularly useful for the calculation of the sub-amplitudes, since 
different gauges can be chosen for different sub-sets of gauge invariant diagrams. 

To conclude this Section, we indicate how these color basis generalizes to the case of loop ampli- 
tudes. First of all let us remind that loop amplitudes can be obtained by applying proper dispersion 
relations to tree-level amplitudes, where some of the external particles have been identified and a 
sum over their possible internal quantum numbers performed. We can then obtain the color form 
factors which generalize our construction to loop amplitudes by contracting pairs of color indices in 
the color representations Eqs. (|3.1|) . (|3.9|) . 

As an example, let us consider one-loop corrections to the qqg\ . . . g n process, whose color struc- 
ture is described at tree-level by Eq. (|3.6|) . For simplicity we will take n = 2. At one-loop we can 
have either a gluon contraction (from a qq plus four gluon tree diagram), or a quark contraction 
(from a qqq'q' plus two gluon tree diagram). Let us study the gluon loop first: for this we need 
to consider the color structure of a qq plus four gluon tree-level diagram, (A Q1 A a2 . . . X a4 )n. Up to 
permutations of the indices, we have three possible independent color structures arising from the 
three inequivalent contractions of gluons: 

N 2 - 1 

5 cd {X a X h X c X d ) ri = 5 cd {X a X c X d X% = — — — (X a X b )ii, (3.17) 
MA c A a A fe A% = ^tr(A a A 6 ) - -^(A a A%, (3.18) 

MA a A c A fe A% = ^(X a X%. (3.19) 

Two comments are in order: first, a term of the form <5jjtr(A a A b ), which was absent at tree-level, is 
now generated. It originates from a color configuration in which the color of the quark flows to the 
antiquark through the gluon in the loop, without emitting any radiation, while the two external 
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gluons are emitted from the remaining color line circulating in the loop. Since the graph of the color 
flow is planar and since no trace over internal color lines appears, this term is of order N°. The 
second comment is that each of the sub-amplitudes that correspond to the three color structures 
(and their permutations): 

N(\ a X%, ^(\ a X%, 5 1J tr(X a \ b ) (3.20) 

is gauge invariant. The proof follows the one given above for the tree-level case. Notice that even 
though the first two color structures are proportional, nevertheless they are independently gauge 
invariant. Graphically, they correspond to planar and non-planar diagrams, respectively. 

From the analysis of the diagrams with a quark contraction, finally, we find again the color form 
factors given in Eq. (j3.20|) plus the form factor (A a A fa )jj (from pure quark- loop diagrams). 

In the general case of n external gluons and m external quark pairs the possible color form 
factors can be represented, in a symbolic fashion, by the following expression: 

NP tr(\...\)...tr(\...\) 5 aBl ...S (\...\)^...(\...\) s - s . (3.21) 

If only external gluons are present, the form factors are given by products of traces of A matrices. In 
general the power p is an integer, determined by the degree of non-planarity of the given color flow 
configuration, by the number of closed color lines, and by the number of U(N)-trace subtractions. 
Once again all of the sub-amplitudes relative to a given form factor are gauge invariant. In spite of 
the proliferation of form factors, which form a highly reducible basis for the color space of a given 
process, the possibility of breaking the sum of diagrams into many gauge invariant components 
turns out to be an extremely efficient bookkeeping device to explicitly carry out the calculations of 
complex matrix elements. 
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4 Super symmetry Relations among Amplitudes 



The properties of the color form factors introduced in the previous Section only depend on the 
representation of the gauge group to which the partons belong and to the gauge nature of the 
couplings, while they are not directly related, for example, to the particle's spin. By this we mean 
that the scattering amplitudes for scalar particles transforming as the fundamental representation 
of SU(N), for example, can be expanded into the same color basis - Equation (|3.6|) - as the 
amplitudes for quarks. This expansion will still be gauge invariant and satisfy the important 
properties illustrated in the previous Section. Likewise, amplitudes with fermions transforming 
according to the adjoint representation of the gauge group (as in a supersymmetric Yang-Mills 
theory) can be expanded using the dual basis, Equation (|3.1|) . 

In a supersymmetric theory, in which particles with different spins are related to one another by 
symmetry transformations, the relation between the color structures extends to a relation between 
the sub-amplitudes as well. This proves to be extremely useful in simplifying the calculations for 
multi-particle processes in supersymmetric theories, as different amplitudes are connected by simple 
algebraic identities. In particular, amplitudes with scalars or fermions are much simpler to evaluate 
than amplitudes with gauge bosons, as the number of diagrams and the complexity of the couplings 
are smaller in the first case. 

The general properties of scattering amplitudes in supersymmetric theories were first discussed 
by Grisaru et al. in Reference |HH l¥Tj . The importance of these supersymmetry relations for 
calculations in non-supersymmetric theories was then pointed out in Reference [7H], where it was 
suggested the use on N = 2 supersymmetry for the evaluation of tree-level multi-gluon processes. As 
was noticed in Reference [TS], in fact, the diagrams contributing to multi-gluon processes at tree-level 
are exactly the same in the ordinary Yang-Mills theory as they are in its supersymmetric extension, 
since neither scalar nor fermionic particles are allowed to appear as internal propagators. The 
amplitudes with only gluons can be related through supersymmetry to easier-to-evaluate amplitudes 
with scalar and fermion external states, thus significantly simplifying the calculations. N — 1 
supersymmetry was also employed by Kunszt [SH] for the calculation of six-parton processes in 



In this Section we will illustrate the basic features of this technique, and we will show how to 
efficiently complement it with the color expansion developed earlier. 

Here we will just use N = 1 supersymmetry, rather than N = 2. 

One possible representation of N = 1 supersymmetry contains a massless vector [g ) and a 
massless spin 1/2 Weyl spinor (A ). The ± refers to the two possible helicity states of the vector 
and the spinor. Let Q(r]) be the supersymmetry charge [7j with r\ being the fermionic parameter of 
the transformation. Then Q(i]) acts on the doublet (g, A) as follows^^ SB; 



QCD. 



[Q{ , n),9 ± (p)) 



T r ± (p,r / ) A±, 




(p)] 



T T {p,r]) g 
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I^G?, rf) is a complex function linear in the anticommuting c- number components of 77 and satisfies: 



T + (p,r]) = [T {jp,r))Y = i] tt_(p), (4.3) 

with w_(p) a negative helicity spinor satisfying the massless Dirac equation with momentum p. 
Because of the arbitrariness in choosing the supersymmetry parameter 77, we choose this to be a 
negative helicity spinor obeying the Dirac equation with an arbitrary massless momentum k times 
a Grassmann variable 9. This variable is used to remind us that r ± (p, 77) anti-commutes with the 
fermion creation and annihilation operators and commutes with the bosonic operators. If we use 
the notation introduced in the first Section, we then obtain: 

T + (p, k) = T + (j9, T](k)) = 6(k + \p-) = 6 [kp]. (4.4) 

As a notation, we choose to label the supersymmetry charge Q(rj) with the momentum k charac- 
terising the parameter r\: Q{k) = Q[rj(k)]. 

Because of supersymmetry, the operator Q(k) annihilates the vacuum. It follows that the 
commutator of Q(k) with any string of operators creating or annihilating vectors g ± and spinors 
A 1 * 1 has a vanishing vacuum expectation value. If Zi represent any of these operators, we then obtain 
the following Supersymmetry Ward identity (SWI) 



= ([Q,i[ *i])o = X>i • • • [Q, Zi]--- z n ) , (4.5) 



i=\ i=l 



where (. . .) indicates the vacuum expectation value. If we substitute in equation (|4.5j) the commu- 
tators, we obtain a relation among scattering amplitudes for particles with different spin. General 
features of Yang-Mills interactions, like helicity conservation in the fermion-fermion- vector vertex 
guarantee the vanishing of some of the amplitudes in (|4.5|l . The arbitrariness in choosing the 
reference momentum k for the supersymmetry parameter 77 allows a further simplification of equa- 
tion (|4.5jl . by choosing k to be equal to one of the external momenta. 

As was first pointed out in Reference [HJ| > these relations can be used to prove general properties 
of helicity amplitudes in supersymmetric Yang-Mills theories; at tree level, these properties hold 
for the non-symmetric theory as well. We will here prove some of these properties as an example of 
the use of the supersymmetry relations. In particular, we will show the vanishing of all the helicity 
amplitudes of the kind (± ± . . . ±) and (=p ± . . . ±), where we assume all of the particles as outgoing. 
For the two-to-two scattering processes in Yang-Mills and gravity, these vanishing theorems were 
first proved in Reference |94j . 

Let us start applying the supersymmetry charge to the following string of operators: 
= (\Q,A+g+g+...g+}) = -F~( Pl , k) A(g+,g+, . . . ,g+) + 



+T + (p 2 ,k) A(A+ A+ ...,0+) + ... + V + (p n} k) A(A+, A+). (4.6) 

Since all of the couplings of fermions to vectors are helicity conserving, all of the amplitudes with 
two fermions of the same helicity must vanish, and as a consequence of the identity the first term 
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on the right hand side of Equation (|4.6|) must vanish as well. Therefore maximal helicity violation 
is forbidden in perturbation theory in any supersymmetric gauge theory, and at tree-level in any 
gauge theory. 

To prove the same theorem for the next helicity violating amplitudes, let us consider the following 
identity: 







(\Q, At 92 9t 1) 



-r-(pi,fc) A(g+,g^,...,g+) - r~(p 2 , k) A(A+, Aj, . . . ,g+). (4.7) 

Here we have omitted all of the vanishing amplitudes with both fermions having the same helicity. 
Equation ()4.7|) must be satisfied for any choice of the vector k, and in particular we can then choose 
k = P2, proving that the gluonic amplitude must vanish, or k — p±, thus proving the vanishing of 
the amplitudes with the fermion pair. 

As a first example of non-vanishing amplitudes, let us now consider the helicity amplitude 
(gf , g 2 , gt, • • • , fi 1 ^), with two negative-helicity gluons and n — 2 positive-helicity gluons where all of 
the particles are outgoing. Through the SWI we can relate this amplitude to amplitudes with two 
fermions and n — 2 vectors. Helicity conservation for the fermions implies that only an amplitude 
with one positive- and one negative-helicity gluino can be non- vanishing. In this way equation (J4.5j) 
reduces to: 

r-( Pl ,k)A{A^,g 2 ,Ai,g+,...,g+) + T~{p 2 , k)A{g^, A 2 , A+ g+, . . . , g+) 



- r-( P3 ,k)A(g^,g 2 -,gt,...,g+) = 0. (4.8) 
Choosing, for example, k = p± we therefore obtain the following relation: 

(12) 

A{g^,g 2 ,g+,...,g+) = hjjAfa, A 2 , A+ g+, . . . , g+) (4.9) 

As we said before, the purely gluonic amplitudes for the non- supersymmetric and the supersym- 
metric theory coincide. 

It is very important to notice that the supersymmetry identities and all the relations that can 
be obtained through their use - like the vanishing theorems - hold separately for each of the the 
sub-amplitudes in which we can expand the full amplitude. This can be easily proved on the basis 
of gauge invariance and leading order orthogonality of the color form factors. 

Another important consequence of the expansion in terms of the color structures described in the 
previous Section is the possibility of relating directly amplitudes with a pair of quarks to amplitudes 
with a pair of gluinos. The SWI will then allow us to connect directly amplitudes with only gluons 
to amplitudes with a quark pair. To make explicit this relation, we remind from the previous Section 
that amplitudes with a quark pair (gluino pair) can be written in the following way: 

M n+2 {q,gi,...,g n ,q) = E (\ ai \ a2 ■■■\ an hm q (q, Pu ..., Pn ,q), (4.10) 

{l,...,n} 
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M n+2 {gi, ■ ■ ■ ,g n , A n+ i, A„+ 2 ) = tr (\ ai \ a2 . . . \ a "+ 2 ) uia^, . . . ,p n ,q n+1 ,q n+2 ), 

{l,...,n+2}' 

(4.11) 

where the momenta labeled with q are fermion momenta, and where the sub-amplitudes m q and ttia 
can be found by summing over subsets of Feynman diagrams obtained according to the prescriptions 
introduced in the previous Section. The prime indicates that only non-cyclic permutations have to 
be summed over. 

The main difference between Equation (|4.10|) and Equation (j4.11|) lies in the fact that while the 
quark sub-amplitudes m q always have the two fermions adjacent, there are gluino sub-amplitudes 
with non-adjacent fermions. The gluino sub-amplitudes, furthermore, satisfy the Dual Ward Iden- 
tity, Equation ([3.4)1 . It is easy to prove, by just studying the structure of the relevant Feynman 
diagrams, that the following identity between quark and gluino amplitudes holds: 

m q (q,p 1 ,...,p n ,q) = m A (q n +2,Pi, ■ ■ ■ ,Pn, Qn+i), (4.12) 

where q = q n+2 and q = q n+ \. The complete proof can be found in Reference 

Therefore by just calculating the gluino amplitudes we can automatically obtain the quark 
amplitudes, using the previous identity, and the purely gluonic amplitudes, by using the SWI. Even 
if we are only interested in the quark amplitudes, it may still be nevertheless useful to consider the 
gluino amplitudes with non-adjacent fermions as auxiliary objects, because they satisfy the Dual 
Ward Identity and can help simplifying expressions. 

We will present a complete example of how this works in detail in the next Section, when 
calculating the five parton processes. 
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5 Explicit Results 



In this Section we will illustrate the use of the various techniques introduced up to now with the 
explicit calculation of some four- and five-parton processes in a SU(N) massless gauge theory. These 
results were known since the original papers j2Sl HH ESI EDI E21 UH EH] , but reproducing them here 
will show the simplicity of these techniques and their advantage over the more standard approach. 
At the end of the section we will collect some results concerning 6-parton processes as well, without 
going into the explicit details of the calculations 6 . We hope this Section will be helpful to the reader 
who wants to familiarize himself with the explicit use of these tools. 

5.1 Four Partons 

As was mentioned in the previous Section it is generally convenient to start the calculations from 
processes with a pair of fermions and to use the supersymmetry relations to simply obtain the am- 
plitudes of purely gluonic processes as a by-product. The use of the polarization vectors introduced 
in the first Section and of the dual color basis, however, makes the four gluon calculation so simple 
by itself that it is useful to just start from it. We will then derive the fermionic amplitude by 
using the SWI. In the next subsection, when describing the five-parton processes, we will follow the 
opposite route, as then the fermionic amplitude calculation is considerably simpler than the gluonic 
one. 

We introduce the following notation: q and q are the momenta of quark and antiquark, h q is 
the quark helicity ( the helicity of the antiquark, hq, is fixed by helicity conservation) and i,t are 
the color indices; pi, (i = 1,2,3), hi and a, will be respectively the momenta, helicities and colors 
of the three gluons. All the particles are taken as outgoing, and therefore momentum conservation 
is given by q + q + Y,Pi — 0. 

All of the diagrams contributing to the four gluon amplitude have the following color structure: 



Here we used the normalization conventions introduced in Section 3 . The Feynman rules are given 
in Appendix C. The Feynman diagrams that enter in the calculation of a given sub-amplitude can 
just be found by imposing the condition that they contain the trace of the string of A matrices in 
the proper permutation. For example, when calculating m(l,2, 3,4) the only diagrams which will 
contribute are those drawn in Fig. 6. 

Notice that the first diagram will also contribute to the subamplitudes corresponding to the 
permutations (1243), (2134) and (2143), but remember that for the calculation of the subamplitudes 
only the kinematical part of the Feynman rules has to be used (see Appendix C). There is only one 
diagram with the four-gluon coupling, and that contributes to all 6 the permutations. 

6 For the details on explicit analytical derivations see (441 192| for 4-quark plus 2-gluons, |59l 1811 for 2-quark 
plus 4-gluons and [75] EI3 E31 E3 EH] f° r the 6-gluon processes. For 7-gluons see jSSj using the recursion relations of 
Section 7 . 



abX eXcd 
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Before proceeding, let us classify the possible helicity configurations. As it was shown using 
the super symmetry relations in the previous Section, the helicity amplitudes with all of the gluons 
having the same helicity, and the amplitudes with all of the gluons but one having the same 
helicity are zero. We can prove this independently here by using an explicit representation for 
the polarization vectors of the gluons, as given in the first Section. In fact, assign to the gluons 
with the same helicity the same reference momentum, and in the case of sub- amplitudes of the 
kind m(=F, ±, . . . , ±) fix this reference momentum to be the momentum of the gluon with opposite 
helicity from the others. Then it is easy to see using the identities given in Appendix A that all 
of the products e l ■ e J vanish. Since by the Feynman rules (or dimensional analysis) it follows that 
at tree- level each diagram will contain at least a factor e l ■ e\ it follows that these amplitudes will 
vanish. 

Therefore for four-gluon scattering the only non-zero amplitudes will be of the form ( h+), up 

to permutations of the indices. Let us then consider the sub-amplitude m(l~, 2~, 3 + , 4 + ), with the 
reference momenta for the gluons (1,2,3,4) given by the momenta of gluons (3,3,2,2), respectively. 
For this choice of reference momenta the only non-zero e« • e.j is e\ ■ €4. Therefore the only non-zero 
diagram from Fig. 6 is the first one, which gives explicitly: 

m(l-,2-,3+,4+) = _ 2ig 2 ^J^Pl^ 

_ (12) 2 [34P " (12) 4 (5 2) 

9 s 12 s 23 ~ W (12) (23) (34) (41)' ^ 

where various definitions and properties of the spinor dot products collected in the Appendix - 
together with the kinematical identity Si2=S34 - were used. Notice that even though the diagram 
containing the t-channel exchange vanishes as a consequence of the choice of reference momenta 
(i.e. the gauge choice), the t-pole (I/S23) appears from the normalization of the polarization vectors, 
signalling that the gauge chosen is singular for s 2 3 = 0. This should be expected since, as we pointed 
out in Section 2 , these gauges are light-cone gauges. Needless to say the result in Eq. (|5.2jl is gauge 
invariant. 

By using the Dual Ward Identity and the invariance under cyclic permutations we obtain the 

following identity, which allows us to express the other inequivalent sub-amplitude m(-\ 1 — ) in 

terms of the one we just evaluated using Feynman diagrams: 

m(l~2 + 3~4 + ) = -m(r3"2 + 4 + )-m(3^r2 + 4 + ). (5.3) 

Applying the Fierz identity Equation (jA.15|) to this DWI gives finally: 



m(l-,2+,3-,4+) = i# (12)(2 3)(34)(4ir ^ 



4 



which generalizes to 



(IJ) 4 

m(l,2,3,4) = ig 2 - - \ - - , (5.5) 
v ' ' ' ; y (12) (23) (34) (41)' v ; 
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where I and J are the indices of the negative helicity gluons. The full amplitude will therefore be 
given by: 

1 



M(l,2,3,4) = ig 2 (Uf E tr(A ai A a2 A fl3 A 



perm' 



(12) (23) (34) (41) 



ig 2 (IJ) A E [tr(X ai X a2 X a3 X ai ) +tr(X a4 X a3 X a2 X ai )] 1 



(12) (23) (34) (41)' 

(5.6) 

where the prime indicates sum over non-cyclic permutations, and the double prime indicates sum 
over permutations up to cyclic and reverse (i.e. (1, 2, .., n) — > (n, .., 2, 1) ) re-orderings. 

In squaring the four gluon amplitude and summing over colors the 0(N~ 2 ) terms in equa- 
tion ()3.5j) can be shown to vanish by using only the general properties, especially the Ward Identity, 
of the sub-amplitude (see Appendix D for the details). Therefore, 

E l-^4| 2 = N 2 (N 2 -l) £ |m(l,2,3,4)| 2 , (5.7) 

colors perm' 

and the square of each sub-amplitude is very simple because the spinor product is the square root 
of twice the dot product. The final result is the standard four gluon matrix element squared: 

E E = N 2 (N 2 - 1) g* (E 4-) E - s \ - • (5-8) 

hel. colors \*>J / perm' 12 23 34 41 

Here we have not averaged over helicities or colors. 

To obtain the amplitude for two gluons plus a qq pair we can use the SWI explicitly given in 
the previous Section for amplitudes of the kind M(=p =p ± . . . ±), Eq. (j4.9|) . and we simply get: 

m{m9 ^ = 19 w¥mwr (5 - 9) 

where we chose the quark and gluon / (7=1,2) to have negative helicity, and where (ql) is a 
short-hand for (qpi). For a negative helicity anti-quark (i.e. positive helicity quark) it is sufficient 
to exchange q with q in the numerator. The full amplitude will be: 

M(q gi g 2 q) = ig 2 (qI) 3 (qI)J2^ ai ^ 2 h r- w iwiowoi ( 5 ' 10 ) 

{ T^ } (<?<?> (<?1> (12) (2g) 

with the following square, summed over colors and helicities: 

1 



E E \ma 2 = 2N(N 2 -i) 9 * e **»% + *%**] E - 

hel. colors \i=l,2 / {1,2} S qq S 12 s ql S (?2 

+ 0(l/N 2 ). (5.11) 
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For the details of the squaring and the explicit form of the sub-leading piece, see the Appendix. It 
is straightforward to check that these results agree with the standard calculations. 

5.2 Five Partons 

We will begin from the calculation of the matrix elements for the scattering of one qq pair and three 
gluons. First of all we classify the possible color form factors for the process. According to Eq. (j3.9j) 
these are given by the 6 permutations of the expression (A ai A a2 A a3 )jj. To these six permutations 
there will correspond six (a priori different) sub-amplitudes. We will consider now the permutation 
(1, 2, 3) of the gluon indices, and will show afterwards how to obtain the others by using the various 
identities introduced previously. Having chosen a color factor, we need to find all of the diagrams 
which contain this given color factor. For the process under study, these diagrams are shown in 
Fig. 7). Notice that only the diagram (a) contributes exclusively to this sub-amplitude. In fact it is 
easy to see that, for example, diagram (b) will also contribute to the sub-amplitude corresponding 
to the permutation (2,1,3), diagram (d) to (1,3,2), (2,3,1) and (3,2,1), and diagram (/) will 
contribute to all six the permutations. According to our technique, in the calculation of a given 
sub- amplitude we will just sum up the terms of each diagram proportional to the corresponding 
permutation. 

Now we have to classify the various possible helicity configurations. Up to permutations and 
charge conjugation, we have four different cases: we can have either all of the three gluons with the 
same helicity as the quark, or just two, or one, or none. As in the purely gluonic case, amplitudes of 
the type M(=p ± . . . ±) vanish identically, as was proven using supersymmetry. We will now prove 
this by using an explicit representation for the polarization vectors. 

Let us consider the case where all the gluons have the same helicity, opposite to the helicity 
of the quark. Let us choose the reference momentum of the gluon polarizations to be the quark 
momentum. It then follows from Eqs. (jA.18|) and (|A.23|) : 

(q±\F{ Pi ,q) = 0, e l -e j = 0. (5.12) 

The bra spinor represents an outgoing quark with helicity ±. Let us then study the branch of gluons 
starting from the first gluon emitted by the quark leg. The only vector quantities that can contract 
with the 7 matrix present at the vertex are the polarization vector of one of the external gluons 
emitted by this branch, or some combination of momenta of the external gluons themselves. In the 
first case the diagram is zero because of the first equation above. In the second case, possible only if 
the branch has more than one external gluon, the saturation of the indices and the dimensionality of 
the couplings (i.e. there can only be at most one power of momentum for each gluon vertex) forces 
at least one scalar product between two polarization vectors. In this case the diagram vanishes 
because of the second identity above. This proof of course extends to tree-level processes with a 
qq pair and an arbitrary number of gluons, and can be easily repeated for the case with all of the 
gluons having the same helicity, equal to the quark's one. 

Let us now consider the case with two gluons (say 1 and 2) having the helicity opposite to 
the quark (say +). According to the matching rule discussed in Appendix B, we will choose the 
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reference momentum for the polarizations of gluons 1 and 2 to be the quark momentum q, and the 
reference momentum for the gluon 3 to be the antiquark momentum q. We then have the following 
identities: 

(q-\t + (Pi,2,q) = 0, e~{p 3 ,q)\q-} = 0, e 1 ■ e 2 = 0. (5.13) 

Using these identities it is straightforward to show that the only non vanishing diagrams are (d) 
and (e) . The evaluation of these two diagrams is very simple if use is made of the various identities 
given in the Appendix A, and leads to the following result for the sub-amplitude: 

m{<r,9T,9t,9t,P) = ^ < g g>< g lt<U g <^><3g> - (5 - 14) 

Before giving the expressions for the other permutations and helicity combinations, we will use 
Eq. (j5.14|) and the supersymmetry transformation to derive the sub-amplitudes for the five gluon 
process. The supersymmetry relation that we need is the following: 

r~(pi,fc)m(Ar,</2>A3>04>05") + r~(p 2 , k)m(gi,A^,At,gt,gt) 

- r~(p 3 ,k)m(gi, g2,g£,gt,g£) = o. (5.15) 

with r~(p, k) = (pk) and A being the fermion field. The second sub-amplitude entering this identity 
corresponds to the quark amplitude we just calculated, as was mentioned in the previous Section, 
while the first term is one of the fermionic sub-amplitudes that would be necessary for the calculation 
of amplitudes with fermions in the adjoint representation. By choosing k = p± we can exclude this 
term, and using Eq. (j5.14j) we directly obtain: 

(12) 4 

m(gl,g;,gt,gt,gt) = igl {12){23mm{51 y ( 5 - 16 ) 

To get the sub-amplitudes for the other permutations, we only need to use the symmetry of the 
sub-amplitude under exchange of identical bosons and the repeated application of the Fierz relation 
(Eq.flSHIIlO on the Dual Ward Identity : 

m(l,2,3,4,5) = 0, (5.17) 

q/ci(l,..,4) 

the sum being over the 4 cyclic permutations of (1, 2, 3, 4). One easily obtains: 

92, gs, <74, g 5 ) = * g s {12){23){M){45){51 y ( 5 - 18 ) 

where / and J are the momenta of the two gluons with the same (negative) helicity. By using 
again now Eq. (j5.15|) one then obtains the expression for the general permutation of the fermionic 
sub- amplitude: 

™(<i-,gug2,gs,t) = ^l ^^^^y (5.19) 
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where now I is the index of the only gluon with negative helicity. Similarly, the sub-amplitude for 
the helicity configuration with one negative helicity gluon and a negative helicity antiquark is given 
by: 

m(q + ,9i,92,9s,r) = ^ ^^^^j - (5-20) 

All of the sub-amplitudes for the processes with opposite helicities (i.e. (+ H )) can be 

obtained from the previous expressions by replacing ( ) products with [ ] products. 

Squaring the full amplitude and summing over colors and helicity configurations, we then obtain: 

\M( 9l ,...,g 5 )\* = 2^Ar3 (iv2 _i)£4 £ 1 > (5 . 21) 

j>j S12S23S34S45S51 

\M(q,q, gi ,g 2 ,g 3 )\ 2 = 2«? S 6 7V 2 (7V 2 - 1) 5>* a * + 8qi 8%) 



]T + 0(N- 2 ). (5.22) 

{1,2,3} s qq s qi s iz s 23S3q 

For the details of the squaring of the color part, see Appendix D. 
5.3 Six Partons 

The six-parton processes are more complex: two independent sets of helicity amplitudes are needed: 
M 2 _4+ and -^3-3+. The first ones are a trivial generalization of the five-parton amplitudes, and 
are given in the case of two quark-four gluon and six gluons, respectively, by: 

M(qt,q;,g;, 9 :,...,g+)=ig 4 (23) 3 (13) £ (A 3 A 4 A 5 A 6 ) 2i 1 (5.23) 

{3,4,5,6} \ 1Z )\ Z6 ) • • • \ Di / 

M(gt,g;,g£,gt,...,gt)=i9 A W E tr(X ai X a2 ... A" 6 ) 1 , (5.24) 

{1,2,3,4,5,6}' \ 1Z )\ Z6 ) • • • \ Di / 

These sub-amplitudes can be shown to satisfy all of the required properties, such as the SWI, 
the Dual Ward Identity and the proper soft and collinear factorization (see Section 6 ). At the 
leading order in iV the sum of these matrix elements squared, summed over colors and over the 

configurations with helicities ( h + ++) and (+ H ), can be easily obtained using the 

properties of the A matrices, giving: 

|M(^,...,^ 6 )| 2 = 2^ 4 (iV 2 -l)£4 £ + 0(iV- 2 ), (5.25) 

t^j s 12 s 23 • • • s 61 

\M(q,q,g u g 2 ,g3,g4)\ 2 = 2g 8 s N 3 (N 2 - 1) Ei(4 s * + s qt s%) 



Eu 2 3 41 1 + 0(N~ 2 ). (5.26) 



SqqS|jlSl2S23 s 34S4g 
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Notice that contrarily to the 4- and 5-gluon case, here the 6-gluon amplitude squared has a non- 
vanishing contribution at the sub-leading order in N. Its precise form is given in Appendix D. 
Using the factorization properties of the amplitude, however, it is easy to check that this sub- 
leading terms do not have collinear divergencies The absence of these enhancement factors 
makes the numerical value of these sub-leading terms even smaller than what one would naively 
expect from the simple 1/N 2 suppression. This fact will be discussed in more detail in Section 10 , 
where we will illustrate some techniques to approximate the multi-parton matrix elements. 

The six-parton helicity amplitudes .M3-3+ is described by three distinct sub- amplitudes, char- 
acterised by three inequivalent helicity orderings: (+ + H ), (+ H 1 ) and (H 1 1 — ). 

Because of duality, as explained in Ref . [oT)] . all of these sub-amplitudes can be written in the fol- 
lowing form: 



m(l,2,...,6) = %g A 



Pi P2 



^123"Sl2-523'545'S56 £234S23"S34<S 56-361 



£345 S 34S45 S 61 S 12 



+ 



■5l2'523 s 34'S45' S 56' S 61 



(5.27) 



where t ijk = {p { + Pj + p k ) 



Sy + Sjk + Ski- The coefficients Pj will depend on the particular 



helicity configuration and on the process (6-gluons or 2-quark plus four-gluons) . For the purely 
gluonic case, a further relation can be found between the P's that will reduce Eq. (j5.27J) to pj] : 



W3+3- (gi,92, ■■■,96) 



a 



P 2 



+ 



r 



^123 s 12'523' s 45 s 56 ^234-S23 s 34'S56 s 61 

£123/^7 + ^2347« + h i5 aP 



t345-534S45S6lSi2 



+ 



(5.28) 



For reference, we give the coefficients Pj's and a, /3, 7 in Table IS~31l5.3l and Table I5~3l respectively, 
without derivation. Here we will just show how to relate the two sets of coefficients, for the purely 
gluonic and the qq plus gluons case, using the various identities introduced in the previous Sections. 

For simplicity we will just work with the ( h ++) helicity ordering, but the same construction 

can be repeated for the other orderings as well. 

Suppose we have calculated the fermionic amplitudes; then it is easy to prove the following 
identity, using a proper SWI: 



[36]m(0! ,g 2 ,g 3 ,gt,9s>9e) = -[ 31 M A 6"j a i ,92 >9s i9t,gt 



-[32}m(Ai,g^,A^,g^gi,g+). (5.29) 

Here by A we refer to a generic fermion, q or q. Helicity conservation has been used to cancel the 
two amplitudes with two negative-helicity fermions, and the Grassmannian nature of T^ 1 was used 
when moving it through A2. The amplitude with the non-adjacent fermions can be extracted by 
using the Dual Ward Identity obtained by moving the gluon 1 : 

m( A i,92, A 3,g4,gt,g£) = 
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-m(A+ , A 3 , g 2 , g 4 , g+ , g+) - m(A+, A 3 , g A , g 2 , g+, g+) 

-m(A+, A3 ,gl ,g+,g 2 ,g+) - m(A+ A3 , g^ , g+, g+, g 2 ). (5.30) 

Therefore the knowledge of the fermionic amplitudes is completely sufficient to obtain the purely 
gluonic ones without having to calculate any additional Feynman diagram. In particular, if one 
were just interested in the numerical value of the amplitudes to calculate scattering processes, one 
could just use the previous equations as operative definitions of the gluonic amplitudes, without 
having to go through the algebra necessary to find explicit expressions. 







Pi 


P2 


P 3 


(93 


#4, #5, #6) 


U =pi + P2 + P3 


V = p 2 + P3 + Pa 


W = p 3 + p A + p 5 


(- 


-+,+)(/) 


[56] 2 <13><23><l|f/|4> 2 





-[16][26](34> 2 (5|^|2> 2 


(+, 


+ ,-,-)(//) 


-[13] [23](56> 2 <4|C/|2> 2 





[34] 2 (16)(26)(1|W|5) 2 


(- 


+ ,+,-)(///) 


[45] 2 <13><23><l|f/|6> 2 


[15](23)(5|\/|3)(4|\/|6) 2 


[45] 2 (16>(26>(1|W|3> 2 


(+, 


-,-,+) (IV) 


-(45> 2 [13][23](6|f/|2> 2 


[16](24)(6|^|4)(3|^|5) 2 


-(45) 2 [16][26](3|1^|2) 2 


(- 


+ ,~,+,){V) 


[46] 2 (13)(23)(l|f/|5) 2 


[16](23)(6|^|3)(4|y|5) 2 


-[16][26](35> 2 (4|^|2> 2 


(+, 


-,+,-,)(vi) 


-[13][23](46) 2 (5|f/|2) 2 


[15](24)(5|\/|4)(3|\/|6) 2 


[35] 2 (16)(26)(1|W|4) 2 



Table 2: The numerator functions Pj for m(gj t ", q 2 , g 3 , g 4 , g 5 , g 6 ). The left column contains the 
helicity orderings of the gluons we define (I\K\J) = (I + \K ■ 7I J+). 
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{93,94,95,96) 


Ps 


(-, -,+,+)(!) 


(23) (34) [56] [61] ((l|£/|4)(2|y|l)(5|W|2) - S 56 [23](34)(5|W|2) - S 34 [56](61)(l|£/|4)) 




[23][34](56)(61)(4|[/|2)(l|y|2)(l|W|5) 


(-)+)+)-)(///) 


-t 123 [15][45](13)(26)(4|y|6)(l|^|3) -t 234 [45] 2 (13)(26)(l|?7|6)(l|^|3) 
+t 3 45[15][45](13)(23)(l|C/|6)(4|F|6) + [45][56](12)(36)(l|[/|6)(4|y|6)(l|W|3) 


(+,-■>-■>+) (IV) 


ti 23 [16][26](24)(45)(3|l/|5)(3|VF|2)+t 23 4[13][26](45) 2 (6|?7|2)(3|W|2) 
-t 345 [13][26](24)(45)(6|C/|2)(3|F|5) + [12][36](34)(45)(6|C/|2)(3|y|5)(3|W|2) 


(~,+,-,+)(V) 


-ti23[16](35)(6|F|3)(4|F|5)(4|W|2)-t 234 [46](35)(6|F|3)(l|C/|5)(4|W|2) 
+t 345 [46](23)(6|y|3)(l|C/|5)(4|y|5)-[46][56](34)(35)(l|C/|5)(4|F|5)(4|W|2) 


(+,-,+,-)(vi) 


[12] [23] [15] [35] (14) (24) <26) (56) <5|Z7|2) + (S 12 S 23 - Si 2 S 45 )[15][35](24)(46)(3|F|6) 
-S 23< S 16 [35] 2 (26)(46)(1|^|4) - [15] 2 [23][34](12)(16)(24)(46)(1|^|4) 
+5 23 5i 6 [15][35](46) 2 (3|l + 5|2) + S 12 S 15 [15][23][35](24)(26)(46) 



Table 3: The numerator functions P s for m(qi, q 2 , g 3 , g 4 , (fe, Qq) with the same notation as Table 1. 



The squaring of the amplitudes is independent of the particular helicity configurations, and the 
explicit formulas for leading and sub- leading terms are given in Appendix D. The same consideration 

concerning the collinear finiteness of the sub-leading piece made before for the helicities ( h+++) 

also holds here. 
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1+2+3 + 4~ c >~fi~ 
i- — ' .) \j 


1 + 2+3~4 +c >~fi~ 

_L z_ t: \j 


1 +2~3+4^ c ) + fi _ 

X £j O T: *J \J 




X = Tit -A- Ttn -\- T)o 


V = T)t -J- T)o -1- T)a 


Si = 71-1 A- T)o -\- Ttr: 


(X 


o 


L J \ / \ 1 \ tJ / 


[x<_)j \<J\u 


(3 


[23](56>(1|X|4> 


[24](56>(l|y|3> 


[51](24}(3|Z|6) 


7 


[12](45>(3|X|6> 


[12](35><4|Y|6) 


[35]<62><1|Z|4> 



Table 4: Coefficients for the ^3+3-(#i, #2, #3, #4, #5, .Ste) sub-amplitudes. We define (I\K\J) = (I + 
\K-j\J+) 



6 Factorization Properties of Dual Amplitudes 

One of the most important properties of the dual amplitudes, which partly accounts for the relative 
simplicity of their explicit expressions, is their factorizability on multi-particle poles. The residues 
at these poles are determined by unitarity, and can be expressed in terms of dual amplitudes for 
processes with a smaller number of external particles. The possibility of factorizing these ampli- 
tudes into products of amplitudes and near-the-pole propagators, puts such severe constraints on 
the amplitudes themselves that often it is possible to deduce their explicit form by just imposing 
unitarity and Lorentz invariance. Subtle cancellations which usually are made explicit only at the 
matrix element square level for the full amplitude, here are made manifest at the matrix element 
level for each single dual amplitude. From the technical point of view, the constraints imposed by 
factorizability provide furthermore a powerful check all along the way while performing complex 
calculations. 

A very simple and instructive way to prove these factorization properties [B7 is by using the 
Koba-Nielsen representation for the amplitudes [531 ESI- While this representation may not be 
too helpful in carrying out explicit calculations 7 , this compact symbolic representation provides 
a powerful tool for deriving general properties of the amplitudes. It was used independently by 
Lipatov in Ref. to study the emission of soft gluons and gravitons, in Ref. jHS] to study 

the production of gluons in tachyon-tachyon scattering and by Fadin and Lipatov in Ref. [S] to 
describe multi-gluon production in a quasi-multi-Regge kinematics, in which all the pairs of final 
state paticles except one have large invariant mass and fixed transverse momentum. 

The following factorization properties can also be proved in a simple and effective way ^2] 
by using the recursive relations introduced by Berends and Giele jTTj and reviewed here in the 

7 The calculation of the five gluon amplitudes has however been carried out explicitly using the Koba-Nielsen 
representation [23 EE]- 
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following Section (see also Ref . |57| for a derivation of the recursive relations using the Koba-Nielsen 
representation of the amplitudes, and further applications of this approach). 

For the sake of definiteness we will deal in this Section with gluonic amplitudes only. As was 
mentioned previously, an n-gluon dual amplitude can be represented by considering the terms with 
lowest momentum dimensionality ([P] 4- ™) in the expansion of the following expression 8 : 

„ n— 1 

m„(&i, e 1 ;k 2 ,e 2 ;...;k n ,e n ) = / JJ dZi /i K n JJ (Zi - Zj) k%kj 

Jz 1 <z 2 <...<z n i=3 n>i>j>l 



expjv- 7 €l€j - + - kl€] V (6.1) 

where \ikn = (z% — z\){z n — z\){z n — z 2 ) is the measure that makes the integral invariant under 
Moebius transformations. The values of z\^ n can be chosen arbitrarily, but usually as follows: 
Z\ — 0, z 2 — 1 and z n = oo. The gluon amplitude is given by the terms in the expansion of the 
Koba-Nielsen expression which are multi-linear in the polarization vectors 6j. 

The singularities of the matrix elements arise from the regions of integration where two or more 
z's coalesce. This follows easily from Eq. (|6.1|) : for example, it is easy to check that poles like 
1/ {hi + kj + . . . + ki) 2 arise from the region of integration Zi ~ Zj ~ . . . ~ z\. From this it follows 
that for a given dual amplitude, represented by a determined permutation of the indices, the only 
singularities that can appear are multi-particle poles in which the indices of the momenta have to 
appear consequently within the given permutation. 



6.1 Soft Gluon Factorization 

Let us start from the simplest kind of singularities, i.e. those due to the emission of soft gluons. 
We want to show that when one of the gluons becomes soft {i.e. p — > 0) the dual amplitude can be 
written as the product of a dual amplitude describing the process involving the remaining gluons 
times an overall factor. 

Let us introduce the following conventions: we will indicate with w the z coordinate of the soft 
gluon, with p its momentum and with ( its polarization. We will take the permutation in which 
the soft gluon is, by convention, inserted between gluon 1 and gluon 2. We will furthermore fix the 
values of zi^, n as given above, and therefore will integrate the soft gluon 'coordinate' w in the range 
Z\ — < w < z 2 = 1. 

It is then simple to prove that in the p — * approximation the Koba-Nielsen formula becomes: 
m n+1 {ki,e 1 ;p,C;k 2 ,e 2 ;...;k n ,e n ) = f dw w pkl (1 - w) pk2 f ZE 

JO J z 1 <w<Z2<-..<z n 

8 For simplicity in this Section we will omit the coupling constant. 
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exp 



eiC , e 2 C hC k 2 ( 



+ 



where 



w 2 (1 — w) 2 w (1 — w) 



n-l 



^ = n n ~ 1 



i=3 



n>i>j>l 



(6.2) 



(6.3) 



ex p Eo 



v£j ^ (^ ^i) 2 



+ 



ki€j 



(6.4) 



The momentum p was kept only in those terms which can give rise to singularities. By expanding 
at the linear level in the polarizations, we will now find integrals in w of the following form: 



I(a,b) 



dw w pkl (l-w) pk2 w- a (l-wY 



T(-q + 1 + pfci)r(-6 + 1 + pk 2 ) 
T(pk\ + pk 2 + 2 — a — b) 



(6.5) 



where the pair a, b can take the following values: (a = 0, 1, 2; b = 0) or (a = 0; b = 0, 1, 2). The 
only integrals which give the leading infrared singularities are 1(1, 0) and 1(0, 1), which in the soft 
limit behave, respectively, like l/pk\ and \jpk 2 . Therefore the dual amplitude corresponding to the 
emission of a soft gluon takes the following form: 



m n+ i(ki, ei;p, C; k 2 , e 2 ; . . . ; k n , e, 



pk 2 pki 



77^(^1, ei;fc2,e 2 ; 



) k n , 6 n ) 



C ' jeik Tft>ni 



(6.6) 



where j ei k is the classical gauge invariant eikonal current. Because of gauge invariance, we can use 
the spinorial representation of the polarization ( with an arbitrary reference momentum - say k 2 . 
For a positive- helicity soft gluon, we find the following result: 



C-3 



eik 



(p\ h [fci) 
V2(k lP )(pk 2 ) 

which is the square root of the usual eikonal factor: 



V2 



(12) 



(lp)( P 2) 



(6.7) 



IC-j 



eik | 



(12) 



(pl)(p2) 



(6.8) 



For the emission of a negative-helicity gluon we just have to change the (ij) products with [ji] 
products. The factorization of the sub-amplitude, Eq. (j6.6|) . does not imply the eikonalization of the 



34 



full matrix element, as is the case in QED, because of the convolution with the color Chan-Paton 
factors: the interference of gluons in the non-abelian theory persists in the soft limit. Repeated 
applications of Eq. (j6.6J) lead to the multi-gluon amplitudes introduced in [§]. The properties of the 
sub-amplitudes in presence of soft-gluon emission were also studied in detail by Berends and Giele 
in Ref . ^2] • Here expressions were given for the case of multiple soft emission, in the case were soft 
gluons are strongly ordered in energy (Ex » E 2 » . . .), and in the case in which energies are not 
strongly ordered (E x ~ E 2 ■ ■ ■ ~ E k « E k+ i . . .). We refer to that paper for the details. 

6.2 Factorization of Collinear Poles 

In a similar fashion one can analyze the factorization properties of the amplitudes near a collinear 
singularity by studying the residues of the appropriate poles in the Koba-Nielsen variables. To this 
end, we will assume that the collinear pair is formed by the first two gluons, and will label the 
variables in the following fashion: the first two gluons will have momenta p\ and p 2 , respectively, 
and polarizations £1,2- Their Koba-Nielsen variables will be w± and w 2 . For the remaining n gluons 
we will use the notation /c i; e« and for momentum, polarization and Koba-Nielsen coordinate, 
respectively. Furthermore we will fix the range of the Koba-Nielsen integration as follows: 

w;i = < w 2 = w < zi = 1 < ... < z n — > 00. (6.9) 

The collinear singularity - l/(pip 2 ) - will arise from the region w — * 0. To isolate the leading 
contributions, therefore, we will expand the KN integral in a Laurent series in w, keeping only the 
singular part: 

m n+2 (p u p 2 ,k l7 ... 7 k n ) = Jo'dW^inti^ (l-^E 2 ^) exp(-PEf) 

! i=2>n _ x [(C1C2K + (PiOOCf - (paCOtf] £(f + f )} , (6-10) 

where E and Z were defined above, and where P = p\ + p 2 . We left out terms like (P1C2XP2C1) 
because they have higher dimension (i.e. they would disappear in the zero slope limit, in the string 
theory language). The integrals in w can be regularized by introducing a factor (1 — w) e , which 
allows them to be defined in terms of Euler functions by analytic continuation, and then taking the 
e — > limit. In this way only the integral in u^P 1 ? 2-1 ) contributes to the leading behaviour. 

By performing the integrations and keeping only the leading terms, we obtain the following 
expression: 

m n+2 ( Pl ,p 2 ,k u ...,k n ) = fl\ i= 2,n-idziZE n™ = i^f fcl ^k) exp(-PEf) 

{(Cic 2 )p,e(S + ¥) + [(CiOQ" + 2(pic 2 )cr - 2(p 2 Ci)c 2 i ni + ¥)} , (6.n) 

where Q = p 2 — p\. The term proportional to (CiQzjPn corresponds to the coupling of a gluon 
with polarization proportional to its momentum. By gauge invariance, after we integrate over the 
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remaining Koba-Nielsen variables it will be proportional to P 2 , and will only contribute to finite 
terms, so in the leading pole approximation we can drop it. What is left can be written in the 
following way: 

1 d 

m n+2 (pi,P2,h,---,k n ) = — rV lx —m n+ i(P,k u ...,k n ), (6.12) 

2(pip 2 ) <9C M 

where : 

V, = [(GC 2 )Q M + 2(pxC2)Cr - 2(p 2 Ci)C 2 1 (6.13) 

is the usual three-gluon vertex, and ( is an 'auxiliary' polarization assigned to the gluon of momen- 
tum P. 

If we select an explicit representation for the helicities, and reintroduce the coupling constant 
(using the normalization conventions given in the Appendix) we obtain the following relations: 

-nV2-.3....) -nnP .:>,.., (6.14) 

m(l + ,2',3,...) ^ { ~i 9 f {12) \ — m(P+ 3,...) (6.15) 

lz{l - z) J S 12 

+ C 9[1 - Z? 1121 1 -^(P-,3,...) 



"§12 



;,;(!-. 2-. 3....) 1_ -1T j M^= \ — m(P~ , 3, . . .), (6.16) 

Sl2 



v^ 1 - *) 



where z is the momentum fraction carried by the first gluons. One can easily check that all of 
the subamplitudes given explicitly in the previous Section do satisfy these relations in the collinear 
limit. 

This Equation shows the collinear factorization of the kinematical part of the dual amplitude. As 
for the color part, factorization can be easily verified by noticing that in the collinear approximation: 

m n+2 (pi, p 2 , h, ■ ■ ■ ,k n ) = -m n+2 {p2,Pi,ki,...,k n ) (6.17) 

and that 

tr(AiA 2 . . . A n ) - tr(X 2 Xi . . . \ n ) = if 12c tr(X c . . . X n ), (6.18) 

which is in fact the product of the color factor of the three-gluon vertex times the color factor of 
an (n — l)-gluon dual amplitude. 

The general factorization properties of the gluon subamplitudes are given by 

m(l,2,...,n) P -^° ]Tm(l,2,...,A;,-P A ) m(p-\k + 1, . . . ,n) (6.19) 
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where P = J2i=i Pi- Of course the full amplitude, including the color factor, must factorize. But 
the color factors introduced in section 3 only factorize to leading order in the number of colors, that 
is, 

tr(X 1 X 2 • • • A") = ^tr(A x ■ ■ • X k X ax )tr(X ax X k+l ■ ■ ■ \ n ) + ^-tr(X 1 X 2 ■ ■ ■ X k )tr(X k+1 ■ ■ ■ A"). (6.20) 

However, the 1 /N terms in the full amplitude cancel at the pole because of the Dual Ward Identities 
for the gluon subamplitudes. 

Similar factorization properties also exist for subamplitudes involving quark- ant iquark pairs. 
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7 Recursive Relations 



The color structure for purely gluonic and processes involving gluons and a quark-antiquark pair 
defined in previous sections allows for the reorganization of the perturbation theory in a efficient 
and straight forward manner. The building blocks are color ordered vector and spinorial currents 
defined with a gluon off mass shell, or a quark or antiquark off mass shell, with all other particles 
on mass shell. If you have calculated these building blocks for n on mass shell legs then there are 
recursion relationships, the Berends-Giele recursion relations, ref. jTT], which allow you to simply 
evaluate these currents with (n + 1) on mass shell leg s. This allows for computer evaluation of 
processes with a large number of external particles'^. A detailed and self-contained description 
of the use of recursive relations in the calculation of multi-parton processes can be found in Giele's 
thesis jHZj- 

7.1 Color Ordered Gluon Currents 

From the set of color truncated Feynman diagrams that make up the subamplitude, m(l, 2, . . . , n), 
one can form a color ordered gluonic current by replacing the polarization vector of the n — th 
gluon with the propagator and allowing the momentum of this gluon to be off mass shell but still 
retain momentum conservation. This color ordered gluonic current will be represented by Fig. 8, 
where the dotted line represents the gluon which is off mass shell. This current will be written as 
. . . ,n — 1) and the subamplitude can be reconstructed from this current by multiplying by 
the inverse propagator and contracting with a polarization vector and allowing the momentum of 
this gluon to be on mass shell, 



Of course these currents, J^, are not gauge invariant and do depend on the choice of reference 
momenta chosen for the (n — 1) on mass shell gluons. Also they depend on the helicity of the on 
mass shell gluons. However these color ordered gluonic currents can be used as building blocks for 
gluonic currents with more external on mass shell legs. 

Consider a gluonic current with n on mass shell gluons. Then the off mass shell gluon is attached 
to the rest of the gluons either through a three or a four point color ordered gluon coupling. At 
these vertices the other legs are attached to color ordered gluonic currents with fewer than n on 
mass shell gluons. This can be seen diagrammatically in Fig. 9. Hence, the color ordered gluonic 
current with n on mass shell gluons can be written in terms o f glu onic currents with less than n on 
mass shell gluons. This is the Berends-Giele recursion relation!^ for gluonic color ordered currents 
and algebraically it is written as 



m(l, 2, . . . , n) = {e"(p„) i[P(l, n - l)] 2 J M (1, . . . , n - 1)} |p ( i,„_i) 



= —Pn1 



(7.1) 



where, P(l,n) = £iP»- 



—i 



n-l 



J M (l,...,n) 



P(l,n) 2 



W(P(1, i), P(i + 1, n)) J u (l, . . . , i) J p {i + 1, . . . , 



n) 



i=l 
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n—1 n—2 

+ E E VA puplJ J v (l, ■ ■ ■ ,i) J p (i + 1, . . . , j) J a (j + 1, . . . , n)} (7.2) 

j=i+l i=l 

where the color ordered three and four gluon vertices are, see Appendix C, 

V3" vp (P,Q) = i^ r (g vp {P-QY + 2g pl * Q v - 2g fiV P p ), 
Q 2 

V A pupa = {2g pp g va - g pv g pa - g pa g vp ). (7.3) 

The current with one on mass shell gluon is defined as 

J„(l) = e>0- (7.4) 

The gluonic currents, J^l, . . . ,n), satisfy properties that are similar to the gluon subamplitude, 
m(l,2,...,n). 

1. Dual Ward identity: 

J M (l,2,3,...,n) + J M (2,l,3,...,n)--- + J M (2, 3, . . . , n, 1) =0. (7.5) 

2. Reflectivity: 

J„(l,...,n) = (-1)" +1 J>,...,1) (7.6) 

3. J(i(l, ■ ■ ■ , n) is conserved: 

F(l,n)" J M (l,...,n) = (7.7) 

There are simple analytical expressions for the color ordered gluonic currents if all the helicities 
are the same or if one is different from the others. Of course we must define the reference momentum 
for the gluons. Here the symbol i for the gluons must be expanded to i% where the z-th gluon has 
helicity A and reference light-like momentum k. Then 

M k ' fe ' • • • ' fc j " 9 V2 (kl) (12) • • • (n - In) (nk) 

and 

Mh,2 k ,...,n k ) - ( 1) g y/ 2 [kl][u] . mm[n _ ln][nk y C 7 - 9 ) 

Berends and Giele, ref.^^, give compact expressions for J M (1 =F ,2 ± , . . . , n ) for a given choice 
of reference momenta. Also, Kosower, ref . [SZj , has given a light-cone formulation of these recursion 
relation to derive the sub-amplitudes m(— , —,—,+,•••, +). Kleiss and Kuijf, Ref.jHHj, have used 
these recursion relations to calculate the 7-gluon amplitudes numerically. 



39 



7.2 Color Ordered Quark Currents 



For the subamplitudes involving a quark-antiquark pair and gluons one can define a Quark and 
Antiquark color ordered spinorial current, see Fig. 10, in a way similar to the gluon currents that 
were defined in the last section. We will write the Quark current as U(q, 1, . . . , n) and the Antiquark 
current as V(l, . . . , n, q). The quark-antiquark pair plus gluon subamplitudes can be obtained from 
these currents as follows: 

m(q,l,...,n,q) = (q | + P(l, n))V(l, . . . , n, q) \q +P[hn)= _ q 



= U(q, 1, . . . , n)(-i){q + P(l, n)) \q) \ g +P(i, n )=-~ q (7.10) 

In manner similar to the gluon current, a recursion relation can be written for this color ordered 
Quark current [TT], see Fig. Hand Appendix C, 



— n_1 _ in i 

U(q,l,...,n) = J2U(q,l,...,m)^=YU™ + l,...,n)- ^— — (7.11) 

m =o V2 (q + P(l,n)) 



and for the Anti-quark current 



mi 



V(l, . . . , n, q) = E * -^yMj (l, . . . , TO )y( m + l, . . . , n , q) (7.12) 

=i (g + P(l,n)) V2 



7/7 = 



and where the spinor currents for the zero gluon case are defined to be 

U{q) = u(q), V(q) = v(q) (7.13) 

in Bjorken and Drell notation. 

These color ordered spinor currents can be defined for massive or massless quarks. For massive 
quarks the propagators in the recursion relations Eqs. (|7.HI7.12j) must be modified by adding the 
appropriate mass term. For massless quarks these spinor currents carry a chirality such that 

(1 ± l5 )V(l,...,n,q±) =0, U{q ± A ) ... ) n){l ± T5 ) =0. (7.14) 

Also for the massless case the zero gluon currents are simply 

F(g ± ) = (g±|, V(t) = \qT). (7.15) 

Again there are simple analytic expressions for these color ordered spinor currents when all the 
gluons have the same helicity as the fermion, 



U(a+ 1+ n + ) - - a n {k - 1 ^ + (7 16) 

U{q ,l h ,...,n h ) - 9 {ql){12) ... {nk) , ( 7 - lb ) 

U{q ,L„,...,n k ) - { g) [qim ... [nk] . ( 7 - i7 J 
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and 



If there is one gluon with opposite helicity to that of the fermion, the spinorial currents are 

"gfc](g + l 
[ql][lk] ' 

fc><g-| 
(?1><1*> ' 

\l-)[kq 
[kl\[lq\ 



U(q + , 1^ ) = (7.20) 



- (7-22) 



and _ _ 

Finally, for two gluons with opposite helicity, we have the following spinorial currents, 

- w ^ < fr + 2 ) » < 1+ 'W + i + 2 »- (7 - 24) 

^'• i - 2 ' ) = biiCK+2)' <2+iw+i+2) - (7 ' 25) 

v W ,2 r ^) = (i + uf )i 2 +> irT ^!M__ , 7 . 26) 

and 

m,-,^) = (i + 2 + j)|l + ) (1+ ^g> (2a) . (7.27) 

A straight forward example using these currents is to calculate the sub-amplitude for (g - , 1 + , 2™, q + ) 
process, 

m(g-,l + ,2-,g + ) = U(q~) i(l + 2 + f) V(lt, 27, q + ) 

^ 2 (92) 3 (g2) 



{ql){12){2q)(qqy 



(7.28) 



which is the previously obtained result. The spinorial currents defined here can be used to derive 
many of the results of other section, especially the section involving multiple gauge groups. 
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8 Exact Results for n-Parton Amplitudes 



The interest in exact matrix elements for n-parton processes in QCD was started in 1986 by Parke 
and Taylor [5Uj who realized that certain non-trivial helicity amplitudes in pure Yang-Mills theory 
could be written in a simple closed form. These same helicity amplitudes have been extended to 
processes including quarks, in Ref . [dTI 1513] . and Vector Bosons, in Ref . [TT1 . and are also used as 
the starting point for approximation schemes, in Ref. [HU EH] , for processes involving a large number 
of partons. These Generalized Parke and Taylor amplitudes are the subject of this section. 

The Berends and Giele recursive relations presented in the previous Section are an extremely 
powerful tool to derive expressions for processes with a large number of partons. They have been 
used up to now to obtain the amplitudes for 7- and 8-gluon processes [TSJ ^] and for processes 
with a color-singlet vector boson and up to 5 colored partons |14| . in this case confirming the 
results independently obtained by Hagiwara and Zeppenfeld in Ref. |47j . Unfortunately most of the 
resulting formulae are very complex and hard to interpret. It is very interesting, however, that for 
some special helicity configurations the matrix elements conserve a very simple universal structure 
independently of the number of particles involved. 

In this Section we present a collection of exact results which hold for a specific set of helicity 
amplitudes with an arbitrary number of particles. The interest in these expressions is not just 
academic, as these results can be used as the starting point for the development of very powerful 
approximation techniques that will be described in a following Section. Furthermore they help 
clarifying the dynamical features of hard multi-particle processes and shed additional light on the 
structure of quantum coherence in the radiation of abelian and non-abelian radiation, as we will 
discuss. 

8.1 Helicity Violating Amplitudes 

Let us consider processes with the following helicity structure: 



We will call these helicity structures Maximally Helicity Violating (MHV) 9 . It is straightforward 
to prove that the only singularities the corresponding amplitudes may have are soft and collinear 
poles of the form l/y^ • Pj). Because of factorization, the residues at these poles are fixed, and 
given by the product of a proper Altarelli-Parisi splitting function times an (n — l)-parton MHV 
amplitude: 

-+...,+)- _^_y75) m (-D(- -, + ...,+). (8.2) 
VWi) 

The A-P function that will appear depends on the helicities of the collinear partons. 
As an example, we can take the 5-gluon amplitude, Eq. f)5. 16j) : 



(-,-,+ 



...,+) or (+,+ 



(8.1) 





(8.3) 



Amplitudes with helicity configuration (+, ...,+) and (— , +,...,+) are equal to zero, see Section 4 . 
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If we study the behaviour of this sub- amplitude near the three inequivalent poles S23 — > 0, S45 — > 
and S51 — > 0, we obtain the following factorization relations: 



(5) 1 

m — > - — r- g 



^ /, m/LvLv/ c -,v » ( 8 - 4 ) 



(23) y \ ^(1 - z) y (IP) (P4) (45) (51) 



^ (45) 9 \jz(l-z) ^ (12) (23) (3P) (PI)' (8 " 5) 



1 



m (5) -> — ^-r ■ // 



^ znovliw^m - ( 8 - 6 ) 



(51) * \z{l-z) y (P2)(23)(34)(4P)' 



with : 



i = zP + 0(^j)), j = (l-z)P + 0(^j)), (8.7) 

% and j being the two collinear momenta. As expected, the first terms on the right-hand side of 
Equations ()8.4|) - (j8.6J) are the four-gluon sub-amplitudes, Equation ()5.5j) . and the second terms are 
the square roots of the polarized A-P functions. 

Lorentz invariance and the factorization properties uniquely fix the form of the amplitudes at 
tree level (their squares are just rational functions), and for the MHV amplitudes these constraints 
can be easily solved explicitly with simple 'educated guesses' of what the amplitude might be. More 
formally, these amplitudes can be derived by solving the Berends and Giele recursive relations, which 
for these helicity configurations turn out to be particularly simple |111 

For purely gluonic processes, the MHV amplitudes are given by the obvious generalizations of 
Equations (Q and lETTfll) 10 : 

il%r, g^,gt, ...,£+)= ^9 n - 2 (l2) 4 ^ £ i tr(AxA 2 ■ • ■ \ n ) -^^±— _ (8.8) 

where the sum is taken over the (n — 1)! non-cyclic permutations of the indices. It is easy to show 
that this Equation satisfies all the required factorization properties and the Dual Ward Identities. 

At the leading order in N the square of these gluonic matrix elements, summed over colors, and 
over all the MHV configurations, gives the so called Parke and Taylor Amplitudes [HHj: 

\M( 9l ,...,g n )\ 2 = 2^- 4 iV"- 2 (V 2 -l)£4 £ . (8.9) 

i>j {1,2,. ..,n}' S 12 S 23-5 3 4 • • • S n l 

The overall factor of 2, coming from the sum over ( !- + ■■■+) and (-H . . . — ) configurations, 

is clearly absent for n = 4 

10 In this Section we will only consider the (— , —,+,...,+) helicities; the (+, +,—,...,—) ones can be obtained by 
replacing (ij) products with [ji] products. 
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Recently a universal form was found [57] also for a specific set of gluonic non-MHV sub- 
amplitudes, namely for helicity configurations of the following form and in the indicated order: 

( h + + . . . +)• This form is not as simple as the Parke and Taylor expression, therefore we 

will not display it here and we refer to the original paper for the explicit result. 

By using the Supersymmetry Ward Identity given in Equation f)4.9|) we can now derive the MHV 
amplitudes for processes with a pair of (massless) gluinos or a qq pair: 

Mfor, A 2 ~ A+ gt, ...,<?+) = ^"- 2 (12> 3 (13> £ tr(X 1 X 2 ■ ■ ■ X n ) 1 

{l,2,...,n}' \ LZ /\ Z6 / " ' ' \ nl ) 

(8.10) 



M(q-,g^gt,...,g+,t)=ig n{ -^P- £ (AxA 2 -.-V ' 



We can once again square these matrix elements at the leading order in N, and obtain: 



(8.11) 



\M(g 1 ,A 2 ,A 3 ,g i ,...,g n )\ 2 = 2^«" 4 iV«- 2 (iV 2 -l)E^ 2i 3(4^ + ^ S 



XEll2 nV 1 > (8-12) 

\M(q,g u ...,g n ,q)\ 2 = 2g 2 s n N n ~ 1 (N 2 - 1) £?=i(4 s * + s^fj 

x— En 2 ni 1 • (8-13) 

Sqq ^i i ,A-">™j' S ql S 12 ...S n q \ > 



Even though massless supersymmetric particles do not exist, the m = approximation might 
turn out to be useful if they were discovered to be relatively light on the scale of the future hadronic 
supercolliders, where their properties would be studied in detail. 

Let us now take the amplitude M(Af, A 2 , A3 , g~l , g£, . . . ,g£)- By commuting with the super- 
symmetry operator and properly choosing the reference momentum k we obtain the following SWI: 

(12) 

M s {A+,f4,A3,At,g+,...,g+) = ^M- g (g+,A+,Aa,gZ,g£,...,g+). (8.14) 
By using Equation ()8.10|) we get: 

M 5 (A+,A+,A3,A 4 -,^ 5 + ,...,^)=^- 2 (12)(34) 3 £ tr{X 1 X 2 • • • A n )-— — ±— — -. (8.15) 

perm* \ / \ / \ / 

The MHV amplitude for the scattering of gluons and a pair of massless scalar-quarks is obtained 
from the SWI and the supersymmetry transformations of a chiral superfield^^^ : 

M(0+,0 2 ,^,^ 4 + ,...,^)=^- 2 (23) 2 (13) 2 £ (A 3 A 4 ■ ■ ■ \ n ) l2il — — 1— — . (8.16) 

{3,...,fi} 
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</) ± are the supersymmetry partners of the two helicity states of the quark. The two combinations 
<f) + ±</>~ transform respectively as a scalar and a pseudoscalar under the Lorentz group. For n = 4, 5 
these are the only independent non-vanishing helicity amplitudes for this process. 

It is interesting to notice, even though the authors do not have a clear understanding of the 
deep reasons for this result, that all of these exact matrix elements for the MHV amplitudes with 
gluons, fermions and scalars can be generated in a straightforward way as correlation functions of 
fields of a two-dimensional Wess-Zumino-Witten gauge model with N = 4 supersymmetry [7Hj . 



8.2 Color Coherence 

Let us go back now to the qq plus gluons process, Equation (J8.ll)) : if we put the color factors 
(A ai . . . X an )ij equal to 1 for each permutation of 1 through n, then Equation (j8.11|) gives rise to 
the QED result for the amplitude with one quark-pair and n photons. This can be easily proved 
diagrammatically by observing that diagrams with non-abelian gluon vertices entering the graph 
expansion for the sub-amplitudes cancel in pairs when we perform the sum over permutations. In 
this way the only diagrams left are the QED-type diagrams, with the common trivial abelian color 
structure. This result is independent of the helicity configuration, and for the helicities considered 
above we then obtain: 

M«(h q , Ky) = i (V2e)" {p7}3 jf 7} J2 J nm! s v ^ 

iPPi {i,w*} ip 1 )! 12 } • • • i n Pr 

where 7 is the momentum of the photon with helicity different from the others, /i 7 is its helicity, 
and where the curly brackets stand for the spinor products or their complex conjugates, depending 
on the helicity of the photon: 

{ij} hy= _ = (ij) , {zj} hj=+ = [tj] (8.18) 
The following remarkable identity holds: 

v {pp} = A {pp} , g 1Q1 

Equation 1)8.19)1 can be proved by iteratively using the Fierz identity: 

{pp}{qq} = {pq}{qp} + {pq}{pq}- (8.20) 

Equation (J8.19)) can be thought of as a sort of 'square root' of the eikonal identity. It allows us to 
put equation 1)8.17)1 into the eikonalized form: 

M 7 >„, M = i i^^gpl ft _gL_ . (8. 21 ) 
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Equations (|8.11|) . (|8.17|) and (|8.21|) offer a nice example of the difference between the properties of 
the non-abelian radiation as opposed to the abelian radiation. Let us take, in fact, the square of 
these three expressions, summed over the colors of the quarks and of the gluons (when present): 



VlM n f/i h)\ 2 = g^N^+i iPdfiPa) \ + _L( inter f) 



(8.22) 



E = (v / 2e) 2 "iV (p7)3( f 7) £ 1 + (inter/.), 

(8.23) 

^^(^^l^^e^^^nT^. ( 8 - 24 ) 



where: 



{ij) = 2i-j. (8.25) 

Equations ()8.23|) and ()8.24|) are identical, thanks to the eikonal identity, but we wrote them in 
the two different ways to establish a connection with the expression for the gluon emission. Equa- 
tion ()8.24)l shows that the photon emission is incoherent: the photons only know about their source, 
i.e. the quark line, but they do not know about each other. Up to the overall factor in front, the 
probability for the emission of n photons is just the product of the probabilities for the independent 
emission of each of them. 11 

On the contrary, if we now look at equation ()8.22|) we see that the gluon emission is not incoher- 
ent: gluons know of each other's presence, and the full probability is not a product of probabilities. 
The interference terms coming from the product of different permutations are suppressed by a fac- 
tor of 1/iV 2 ; this suppression originates from the interferences of the color factors. For the photon 
emission, vice versa, we can see from equation ()8.23j) that the interferences among different per- 
mutations are not suppressed and they conspire to cancel the coherence apparent into the sum of 
squares, giving rise to the factorized expression given in (|8.24j) . 

The phenomenological consequences of this coherence effects have been explored experimentally 
jil] (the string effect) and theoretically (see for example El El EH US 123 EH)- 

11 This result, which is exact for this specific helicity configuration, also holds for any other helicity configuration 
in the limit of soft-photon emission. The reason why it cannot hold for an arbitrary helicity configuration is that in 
general the amplitude will have poles of the kind l/(p + k + k 1 ) 2 , k and k' being arbitrary photon momenta. 
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8.3 qqqq plus gluons 

In general the factorization of the color structure exhibited in equation (|3.9|) does not imply a 
similar factorization of the kinematical part of the amplitude. In other words, the sub-amplitude 
that multiplies a given color factor does not factorize into products of terms that only depend 
upon the kinematical variables (helicities and momenta) of the particles belonging to the same 
antenna. One remarkable non-trivial exception to this general feature is given by the amplitude 
for a process with two quark pairs and an arbitrary number of like-helicity gluons (all the particles 
are outgoing). Up to an overall factor that only depends upon the helicity configuration each sub- 
amplitude factorizes into the product of two terms that only depend upon the momenta of the 
gluons emitted by one or the other of the two antennas: 

M(h p ,h q ,h g ) = ig n+2 A {h p ,h q ,hg) ■ 



E 



{pq} 



{qp} 



{pa 1 }{a 1 a 2 } ■ ■ ■ {a k q} {g&i}{6i& 2 } • • • {h'P} 

1 iPP} i^} ( X ai X a k) (xbl X b kl] 

N {pa 1 }{a 1 a 2 }...{a k p}{qb 1 }{b l b 2 }---{h>qr ' W W 



(8.26) 



The helicity structure of this amplitude uniquely determines the pole structure and the residues of 
these poles, through unitarity: Equation (j8.26J) is the only Lorentz invariant amplitude that gives 
rise to the right poles and the right residues. Alternatively, one can prove Equation ()8.26j) by using 
the appropriate form of the Berends and Giele recursive relations jllj . 

If the gluons have all positive helicity, then {ij} = (ij), otherwise {ij} = [ij]. The indices p and 
q (p and q) refer to the quarks (antiquarks) and the indices a a , bp refer to the gluons. The arguments 
h represent the helicities of the two quarks and of the gluons; the helicities of the anti-quarks are 
fixed by helicity conservation along the fermion lines. The sum is over all the partitions of the n 
gluons (k + k' — n, k — 0, 1, . . . , n) and over the permutations of the gluon indices. When k = (0, n) 
the product of zero A matrices becomes a Kronecker delta and one of the two kinematical factors is 
equal to one. The overall factor A can be written as follows: 

The functions clq are given in Table where helicity configurations obtained by permuting the 
quark helicities have been omitted. The functions a are universal, in the sense that they only 
depend upon the spin-1/2 nature of the quarks. As we will see later, they also enter in processes 
like deep inelastic scattering or e + e~ annihilation. 
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(h h h ) 


Lb ()\ IL pi IL qi lb g) 


(+, + , +) 


(pq} 2 [pp\[qq\ 


(+ + ~) 


\t>o]^ (nn) (on) 
ipyj \FF/ \\m/ 


(+,"+) 


(pq) 2 [pp}[qq] 


(+,--) 


[pq] 2 (pp) (qq) 



Table 5: The universal functions ao(h p , h q , h g ). 



To the leading order in N, the amplitude squared summed over colors is furthermore given by: 
-£\M(h p ,h q ,h g )\ 2 = g 2n+4 N n (N 2 -l)\A (h p ,h q ,h g )\ 2 

col 

£ (gg) (gg) / 8 28 ) 

(pai){aia 2 ) ■ ■ ■ {a k q) {qb 1 ){b 1 b 2 ) ■ ■ ■ {hyp)' 

If the quarks are identical we must add the contribution from the crossed channel p <-» q. 

As in the case with one quark pair, we can here compare the properties of photon radiation with 
those of gluon radiation. A reasoning similar to the one used in the previous section allows us to 
write the amplitude for the emission of n like-helicity photons off two quark-pairs: 

M(h p ,h q ,h g ) = ig 2 (V2e) n A (h p ,h q ,h g ) ■ 

^ ww • • • {a k p} sim • • • {m> {5hh5i2ji ~ N 5hh5i2h) - 

(8.29) 



Only the contribution from gluon exchange is shown. The effect of photon exchange between the 
two quark-pairs can be easily added. A repeated use of the Fierz identity, equation ()8.2U|h leads 
then to the following form of equation (|8.29|) : 

M(h p} h q} h g ) = ig 2 (V2e) n A (h p ,h q ,h g ) ■ 

nILl + (8 ' 30) 

This expression shows that photons are emitted independently. Once again we expect this result to 
hold for an arbitrary helicity configuration in the soft-photon limit. 
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If we substitute the color factor in equation (|8.3U|) with the Abelian one, 6^6^, and if we put 
g=^/2e, then we obtain the amplitude for the process e + e~ /i + yT photons, as given in reference [THJ 



8.4 e+e and DIS 

It is easy to derive expressions analogous to ()8.26|) and ()8.28|) for the process llqq gluons, where 11 
is a lepton-antilepton pair (for example e + e~ or e~v) [TT1 164j . Again the gluons have all the same 
helicity 12 : 

M l (h h h q ,h g ) = ig n £ M v (h h h q ,h g ) £ (A ai . . . A ar % TZTTJT^' CT > ( 8 ' 31 ) 

V=~/,Z,W {l,2,...,n} iy i Jl iz J • • • l n( lS 

EI^AA)I 2 = ^^(^-1)1 £ M*y(*fAA)l 2 £ ■ (8.32) 

co( V=~/,Z,W {l,...,n} \" 'V / " " " V "/ 

g and g are the quark momenta and 1 through n are the gluon momenta. The contributions from 
photon, W and Z exchange are explicitly exhibited. The functions M v (hi, h q , h g ) are given by: 



Mv(h h h q , h g ) = ^^g^ a ( hl , h q , h g ) (8.33) 



Qv(hi) {Qv{h q )) is the charge corresponding to the interaction of a lepton (quark) of helicity hi 
{h q ) with the vector V. Furthermore s = (p + p) 2 , with p and p being the lepton momenta, and 
M v is the mass squared of the vector boson V. The universal functions a (hi, h q , h g ) coincide with 
those given in Tabled 

For e + e~ scattering the effect of photon radiation (both from the initial and the final state) can 
be easily incorporated into equation (|8.31|) by using equation (|8.30|) . Here we will display directly 
the result for the square of the amplitude with n gluons and m photons, to the leading order in N: 

]T|M e+e " (h,h q ,h g )\ 2 = (2e 2 ) m g 2n N n -\N 2 - 1) | £ My (hi, h q , h g )\ 2 x 

col V=-y,Z 

TT I I {PP} | to) N|2 V (gg) /oo 4 n 

Li {pki}{hp} + {1 £ in} (gl)(12) . . . W 1 • J 

The fcj's are the momenta of the photons. Once again this result is only exact if all the gluons and 
the photons have the same helicity, but this is the behaviour of all the other helicity configurations 
in the case of soft emission. 

Exact expressions for the full calculation of the processes e + e _ — ► 4 partons were given in [20] 
(for the complete 0(a 2 ) calculation) and in pQ (tree level) and for the processes e + e~ — > 5 partons 
in (mill] (tree level). 



2 In Ref.|M| the term {qq} was inadvertently omitted from these equations. 
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9 Multiple Gauge Groups 



In a theory in which the fermions are coupled to a direct product gauge group, say SU(M) x 
SU(N), then the techniques described earlier in this report can also be implemented. In fact 
the subamplitudes can be easily obtained from the subamplitudes given earlier by adding various 
combinations of 'color' orderings together. This can also be extended to U(l) gauge groups, e.g. 
electromagnetism. For spontaneously broken gauge groups, only processes involving a single vector 
boson are easily obtained from previous calculations. In this section we have set all coupling 
constants to unity, but the reader can easily insert the appropriate couplings. 



9.1 SU(M) x SU(N) 

Consider a fermion which is in the fundamental representation of both SU(M) and SU(N), then 
the amplitude for scattering of a fermion- antifermion with m vectors from SU(M) and n vectors of 
SU(N) can be written as 



M(g,l,2,...,m;l,2,...,n,g) = £ (A 1 A 2 . . . A m )i j (X 1 X 2 . . . X% j 

p,p 

m M ,N(q, 1,2, . . .,m,q; q, 1, 2, . . . , n, q). (9.1) 



A and A are the fundamental matrices of SU(M) and SU(N) respectively. Also ij and ij are the 
SU (M) and SU(N) 'color' indices of the fermions and the sum is over all permutations of 1, 2, . . . , m 
and 1,2, ... ,n. 

The subamplitudes defined by Eq. (j9.1|) can be easily obtained from the subamplitudes obtained 
earlier, Eq. P~TUJ) . 

m M ,N(q, 1,2, ■■■,m,q; q, 1, 2, . . . , n, q) = m (<7, 1, 2, . . . , m, 1, 2, . . . , n, g) (9.2) 

where the J2i is over all ways the barred numbers can be interspersed within the unbarred numbers 
maintaining the order of both the barred and unbarred numbers. This sum causes all Feynman 
diagrams which connect directly the vectors of SU(M) with SU(N) to be cancelled. 

As an example, consider the scattering in which the fermion has negative helicity and the a 
vector boson of either gauge group has negative helicity and all particles have positive helicity. 
Then, 

/ , _ t — _\ . (qa) 3 (qa) (qq) 

m M ,N{q, l,...,m,q;q,l,...,n,q) = i — 2 2^ 



(qq) 2 Y (gl)(12)---(ml)(12)---(fig) 

(qa)\qa) (qq) (qq) , . 

{qq) 2 (ql)(12)---(mq) (ql) (12) • • • (fiq) ' 1 " ' 
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9.2 U(l) 



For completeness we include here a discussion of the U(l) gauge group even though it was extensively 
discussed in the previous section. The U(l) gauge group results can be obtained by replacing A in 
Eq. (j4.10|) with 5 or by iterating Eq. ()9.2j) of the previous section. Thus for QED the amplitudes 
are 

A QE D{q, 1,2,..., m,q) = J2 m ^ ^ 2 , • • • , g). (9.4) 

p 

The sum over all permutations causes all non-abelian Feynman diagrams that naively appear to be 
cancelled, thus leaving only the QED Feynman diagrams. 

Again as an example consider the case of fermion-antifermion m photon scattering in which the 
fermion and the a photon have negative helicity and all other particles have positive helicity. Then, 

A ( , — \ -(qafiq®) (qq) 



(qqy i?(ql)(12)---(mq) 
(qafjqa) n (qq) , , 

~ (m) 2 Y t(Qmy { } 

Combining this example with that of the previous subsection, the subamplitude in an SU(M) x 
SU(N) x £7(1) theory for a fermion and the a vector particle with negative helicity and all other 
particles of positive helicity is 

m M ,N,QED(q, 1, • • • , n, q; q, 1, . . . , n, q; q, 1, . . . , n, q) = 

. (qaf(qa) (qq) (qq) ^ (qq) 



(qq) 2 (ql)(12)---(mq) (ql) (12) • • • (nq) ¥ (qi)(iq) 



IItTO- (9-6) 



9.3 The Insertion of a W or Z 



A spontaneously broken gauge group does not have a simple generalization of the previous subsec- 
tions. However, the insertion of one such massive vector particle, W or Z, can be easily incorporated. 
Consider the scattering of a quark-antiquark n gluons and a W vector boson. Then the amplitude 
for this process is written as 

A(q, 1, . . . , n, q; W) = ^(A 1 • • • X n ) fj m(q, 1, . . . , n, q; W) (9.7) 

p 

where the subamplitude can be written as 

m(q, 1, . . . , n, q; W) = i £ U(q~, 1, . . . , i)l } l ~ lb) V{i + I, . . . , n, g+). (9.8) 

i=0 1 

Here the recursion techniques of section 7 have been explicitly used. 
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This expression can be used in one of two ways: either one can square it directly or allow the 
W boson to decay into another fermion-antifermion pair. If one squares this expression directly the 
relationship 

. , W^W U 

E«W = -9" U + ^7T- (9-9) 

pol 1V1 W 

can be employed. 

The other alternative is to replace the polarization of the W vector boson by the amplitude for 
it to decay into a lepton-antilepton pair 13 . Then Eq. f!9.8|) is written as 

m(q,l,...,n,q;L,L) = -i U{L~) ltl (1 ~ 7s V (I + ) 



X w 



(W 2 — Myy + iM w Y w ) 
x J2U(q~A,...,thv K -—^V(i + l,...,n,q + )- (9.10) 

If we use the fact that the charged lepton is effectively massless, compared to Mw, the Fierz 
rearrangement gives 

m{q,l,...,n,q;L,L) = -2« ^ — — —j— . (9.11) 

(W 2 - Mfo + iM w T w ) 

Using the results from the recursion relation section of this report, the sub-amplitude for the 
process qq — > W — > LL is 

,n w ( r ,f,L-,L + ) - -*F«I<» £ > 



(W 2 - Ml- + iM w T w) 

2i {qL) 2 [LL] 
(qq) (W 2 -M 2 V + iM w T w ) 

2i [qZf (LL) 
[qq] {W 2 -M 2 v + iM w T w y 

Adding n gluons with the same helicity to this process, gives 



(9.12) 



m (a' n+ n+ 77+ • J ~ V \ - 2z (qL) 2 [LL] 



( - -+ T T+^ (-1)" 2i [q L] 2 (LL) 

m w {q ,g 1 ,---,g n ,T)L ,L ) = . 9.14 

[<?1][12] • • • [nq\ (W 2 - Mfr + iM w T w ) 



13 Working at the amplitude level one could also use the representation for the heavy vector polarizations given in 
Ref.[Hl]. 
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If we add two gluons of opposite helicity, then the sub-amplitudes are 



, (qL)(L+\q + W\2+)[lq] 2 
S q w S12 [2q] 



-2i 



+ 



+ 



(W 2 -M 2 V + iM w Y w ) 
(q2)(qL)[Lq][lq] 

(?1> Si 2 [2q] 
(q2) 2 (l+\^ + W\L+)[Lq] 



(9.15) 



(ql) Si 2 Sqw 



and 



m w {q ,9i ,g£,q + ;L ,L 



■+ 



2i 




+ 



+ 



(W 2 -M 2 V + iM w T w ) 
(2+|f +W\L+) (L+\q + W\l+) 

[ql] S12 (2q) 
[q2] (ql) (2+\^+W\L+)[Lq} 

[ql] Su SgW 



(9.16) 



These expressions reproduce the results first obtained in [131 EH The full set of radiative 

corrections (to order a 2 ) to this process was recently calculated in jUEH]- 

The previous discussion can be extended to include the Z boson by decomposing the coupling of 
the Z to the quarks and leptons into its left and right handed parts and then proceeding as with the 
W boson. For a complete discussion of the calculation for these processes, including the complete 
results for processes including a W boson plus five partons see Berends, Giele and Kuijf, ref . |14j . 
These results agree with those independently obtained by Hagiwara and Zeppenfeld, ref. j47j . 
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10 Approximate Matrix Elements 



The techniques described in the previous Sections provide very powerful tools to calculate the matrix 
elements of very complex processes. As an example, the Berends and Giele recursive relations were 
recently used for the calculation of 8-gluon scattering ^B] . The resulting expressions, however, prove 
very slow to evaluate numerically because of their complexity, thus making it almost impossible to 
generate a number of events large enough to perform relevant physics studies. 

These considerations, and the importance of having fast event generators to simulate multi-jet 
processes at high-energy hadron colliders, where these processes will provide important backgrounds 
to many possible new physics signals, justify the study of approximate expressions which describe 
sufficiently well the exact matrix elements throughout phase-space and at the same time are simple 
enough to allow very fast simulations. 

Kunszt and Stirling [UT] and Maxwell [73] were the first to realize that the Parke and Taylor 
amplitudes, Equation (|8.9jl . can be properly fudged in a systematic way so as to reproduce the full 
sum over all the allowed helicity amplitudes for gluonic processes. This idea was later generalized to 
other processes in which at least one set of helicity amplitudes is known in both hadronic [73|Iq^H I74] 
and e + e~ [7S] multi-jet production. An alternative scheme based on the non-abelian version of the 
eikonal approximation was also introduced in in this Section we will describe these various 
approximation schemes, referring the reader to the original literature for numerical comparisons 
between them. 

10.1 The Kunszt and Stirling Approximation 

We will start from the simplest scheme, namely that of Kunszt and Stirling (KS, see Ref.|H° ])■ It 
amounts to assuming that all of the helicity amplitudes have 'on average' the same value, and there- 
fore the full amplitude can just be obtained by multiplying the Parke and Taylor (PT) expressions 
by a proper weight, representing the ratio between the number of non-zero helicity configurations 
and the number of the Maximum Helicity Violating (MHV) configurations whose matrix-elements 
are described by the PT formula. 

This approximation becomes particularly simple when neglecting sub-leading terms in 1/JV. 
This is justified because the sub-leading terms have softer collinear singularities than the leading 
ones, and therefore do not contribute substantially to the numerical value of the matrix elements. 
In particular, for n = 6 the sub-leading terms are finite [HS], and only contribute of the order of 
few percent to the full square. For generic n it was proven in Reference [33] that the sub-leading 
terms are also finite in the strong energy-ordering kinematical domain. This important result 
strongly justifies neglecting the sub-leading terms at the level of precision given by these tree- level 
calculations. 

For an n-gluon process the number of MHV amplitudes is n{n — 1) if n > 4 and n(n — l)/2 if 
n = 4. The total number of non-zero helicity amplitudes is instead 2 n — 2(n + 1). For n = 4, 5 
these multiplicities coincide, and the PT formula describes the exact results, as is well known. For 
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n larger than 5, the KS approximation gives: 



2 n — 2(r> 4- 1 1 

da^ 9 - 9 = n{n [ _^ } daft"-* (10.1) 

For n = 6, 7 , for example, the fudge factor is 5/3 and 8/3, respectively. 

To describe processes with initial state quarks, KS suggest the use of the so called effective 
structure function approximation [2SII1H], which gives a good description of the two-to-two QCD 
processes. According to this approximation in most of the relevant phase-space the differential 
cross-sections for processes initiated by gg, by qg and by qq or qq stand in a constant ratio: 

da gg : da gq : da qq = 1:4/9: (4/9) 2 . (10.2) 

In this way the full differential cross-section, weighted by the appropriate structure functions, reads: 

da tot = F(x 1 )F(x 2 )da gg , (10.3) 
F(x) = g{x)+A/9{q{x) + q(x)), (10.4) 

g(x) and q(x) being the gluon and quark structure functions. For da gg , finally, one takes Eq. (jl(J.l|) . 

The KS approximation scheme tends to overestimate the exact results and the effective structure 
function approximation is works less and less for an increasing number of partons in the final state; 
nevertheless the KS approximation is an extremely useful tool for simple but significant estimates 
of multi-jet rates and distributions. For comparisons of this scheme with exact calculations, see for 
example References jHH EM EH E3 ■ 



10.2 The Infrared Reduction Technique 

It is well known that in the limit in which two partons (say i and j) become collinear, a given 
process can be described in the Weisszaker- Williams (W-W) approximation: 

da™ = — ^ f(z) da^ (10.5) 

where f(z) is an appropriate function of the fraction of momentum carried by one of the two partons 
becoming collinear, and da^ n ~^ is the partonic cross-section for the effective (n— l)-particle process 
in which the two collinear partons are replaced by the single one into which they merge. On the 
pole the W-W approximation is nothing but the factorization of the amplitude, discussed in various 
occasions in the previous Sections. The functions f(z), in the case of a QCD process, are just the 
Altarelli-Parisi (AP) |2. splitting functions. 

The infrared reduction technique introduced by Maxwell [73] improves the W-W approximation 
by using the exact matrix elements for some simple helicity configurations, and derives the other 
helicity configurations by approximating their relative weights at the closest collinear pole. 
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Next to a collinear pole (say pi ■ p 2 — > 0) each of the non-vanishing helicity amplitudes will 
factorize in the following way: 

daP = Mz) dat 1] + finite (10.6) 

where h! are the various helicity configurations which can contribute to the factorization, and fy (z) 
are the corresponding polarized AP splitting functions, depending on the variable z = Ei/(E\ + 
E2). For the time being we will restrict our attention to gluon scattering. For the full process, 
factorization is described by Equation (jl0.5|) . with f(z) given by: 

If we just sum over the PT amplitudes, instead, we obtain: 



dvpT = 7T7 7 fpr(z, s i:j ) da { p T 1] (10. 



where da^ is the sum over all the MHV amplitudes, and fpriz, Sy) is given by: 

f PT (z, Sij ) = 9 2 N R + Z *+} 1 ~ Zr (10.9) 



R = (10.10) 

l^i S Pi 

the indices i and j being different from the collinear particles, and P being the sum of the collinear 
momenta. 

Equation (jlO.Sj) can also be rewritten in the following fashion: 

dof T = x- 1 fAp(z) d*^ (10.11) 



with: 

(l + it!)(l + z 4 + (l -^ 4 
R + z 4 + (1 - z 



x(z, Sij ) = ^ „ A 7 , J (10.12) 



By equating Equations (|10.11|) and (|10.5|) we therefore obtain: 

da [n ~ 1] 

dofuii = da$x{z,Sij) 7! n (10.13) 

d<Tprp 

Maxwell suggested that while the W-W approximation is not very good unless we are very close to 
a collinear pole, Equation (jlO.lHj) is rather good throughout phase-space, provided we perform the 
factorization considering the pair of partons with the minimum s^. In other words, while the value 



56 



of the full differential cross section is not well reproduced in the W-W approximation away from 
the collinear poles, what is well approximated is the relative weight of different helicity amplitudes. 

Since dafl u = dap%, for n = 6 we obtain: 

dvflu = daP TX (z, Sij ) (10.14) 
while for larger n the infrared reduction can be iterated, giving: 

n 

dvfuii = da^Uxkiz^Sij), (10.15) 

fc=6 

with an obvious notation. 

If the two partons which minimize \sij\ belong to initial and final state, we can still use Equa- 
tions f)10.13|) and ()10.12j) provided we keep all of the momenta as outgoing (which implies that the 
energies of the initial state particles will be negative) and define: 

L ' (10.16) 



K + E. 



The z defined in this way cannot be interpreted directly as the fraction of momentum anymore, 
since it will not satisfy the constraint < z < 1. In particular, if i is the final state parton then 
z < 0, while if i is the initial state, then z > 1. However it can be easily checked that with this 
prescription Equations (jl0.13|) and (jl0.12|) reproduce the desired factorization properties. 

This technique can be applied whenever we have exact expressions for some sets of helicity 
amplitudes. In particular, it applies to qqg . . . g processes Elj and to e + e~ — > qqg . . . g and 
eq — > eqg . . . g (75]. We will here summarize the main results concerning the quark-gluon processes. 
Similarly to the purely gluonic case, the infrared reduction technique leads to the following relation: 

n 

da fuU (qqng) = da MHV (qqng) J[ Xk(z k , s^-), (10.17) 

k=A 

where the relevant MHV amplitudes were given in the previous Section, Equation (|8.13|) . For these 
processes the factors \ depend on the nature of the partons that have the minimum |sy|. If these 
are both gluons (with indices a and j3), then we have: 

(1+R) (l+^ + (l-z) 4 ) 

Kaa = ; a ; rr^ (10.18) 

A " (R + z 4 + (1- z) 4 ) V ; 

as before, but with 



P = Pa+Pf3 



R = "ff lViiT 7y . (10.19) 

[SqpSqP + S q pS 6 qp ) 
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If the pair with the minimum dot product contains a quark and a gluon then 



= (10 - 20) 



where 



Q° 



z i=™ > Q = Pa + q 



R = ggglflif|i . (10.21) 

Ei^a S Qi S qi 

The result for an antiquark-gluon pair is the same as the above quark-gluon pair but with each 
fermion momentum replaced by the appropriate anti-fermion momentum. 

For the situation in which the minimum \sij\ pair is made up of a quark and an antiquark the 
multiplication factor is 

Xqi = (l + R) (10.22) 

where 

G = q + q 



R = t^lA, (10.23) 



J2i<j s 



In all of these cases we assume the prescription given above when the collinear partons belong 
to initial and final state. 

The Maxwell approximation scheme has been checked against exact matrix elements for various 
multi-parton processes, and has proved to be extremely accurate, in addition to being numerically 
more efficient by 2 or 3 orders of magnitude, depending on the number of IR reduction steps. 
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11 Conclusions 



While most of the explicit techniques summarized here are limited in their application to tree-level 
processes, it is auspicable that one day they can be extended for use in loop calculations as well. For 
example, the color structures introduced in Section 3 provide a gauge invariant decomposition of the 
amplitude at any order in perturbation theory. The spinor representation for the polarization vectors 
introduced in Section 2 and used throughout this report is however specific to four dimensions and 
could not be used in a dimensional regularization scheme. We believe it can be extended to non- 
integer dimensions in a scheme which preserves the dimensional relations between spinors and 
vectors, such as supersymmetry. However such a scheme (dimensional reduction) may not provide 
a consistent regularization of loop amplitudes beyond one loop. 

Perhaps some of the beautiful features of amplitudes with simple helicity configurations (the 
Parke and Taylor amplitudes) can be shown to persist at higher orders in the loop expansion. 
The structure of multi-parton amplitudes unveiled by the approaches described in this report will 
hopefully lead to a better understanding of perturbation theory for non-abelian gauge theories. 

The intriguing connection between the Parke and Taylor amplitudes and correlation functions 
in a two-dimensional Wess-Zumino-Witten model with N = 4 supersymmetry discovered by Nair 
suggests the existence of new generating functionals for these amplitudes, and perhaps more fun- 
damental structures underlying the perturbative expansion in a gauge theory. 

For purely phenomenological purposes, the production rates obtained by these calculations at 
tree level are extremely valuable for reliable estimates of important processes. These complex 
processes with many partons in the final state, are now being experimentially probed by current 
hadron colliders and will become more important at the next generation of hadron colliders currently 
under construction. A detailed understanding of these QCD/Electro-Weak background processes 
will be fundamental for the detection of signals of new physics which will contain many jets in the 
final state. 
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A Appendix: Polarization Vectors and Spinor Properties 



We will use notations and conventions as in [92], and will summarize them here for ease of reference. 
Let ip(p) be a massless Dirac spinor. We will denote its chiral projections as follows: 

1 



\p±) = V±0) = -(1±7sMp) (p±I = Mp) (A- 1 
By convention, we will choose the spinor phases so as to satisfy the following identities: 

|p±) = \ P T) C (P±\ = %PT\, (A.2 
where the suffix c stands for the charge conjugation operation: 

I > c = C\ >*, c (| = -^l^- 1 (A.3 

Cl^C- 1 = 7 „, (A.4 

C = C ] = C- 1 = C* = C T . (A.5 
We will also introduce the following notation: 

(pq) = (p-\q+) [pq} = (p + \q-} (A.6 
The spinors are normalized as follows: 

(p\lM = 2p„ (A.7 
From the properties of the Dirac algebra, it is straightforward to prove the following useful identities 

(p + \q+) = (p — \q—) = (pp) = [pp] = (A. 8 

ipq) = -(qp), [pq] = -[qp] (A.9 

2\p±)(q±\ = ^(l± l5 )Y(q±\l»\P±), (A.10 

(pq)* = —sign{p ■ q)\pq] = sign(p ■ q)[qp] (A. 11 

\(pq)\ 2 = 2(p-q), (A.12 

(P±|7 W •••7w„ + i |?±) = (9T|7 Wn+ i •••7w \pT), (A. 13 
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(p±|7w-7tfanl9T) = -(9±l7M 2 n---7w IpT), 



(AB)(CD) = (AD)(CB) + (AC)(BD) 



(4 +|7/i (C - 1 7" |I>-) = 2L4D] (C5). 



(A.14) 
(A.15) 
(A.16) 



In the identity Eq. (|A.11|) the possibility of having spinors with energies of different sign is considered. 
This will be important in the following, since for simplicity we will always carry out the calcula- 
tions of the matrix elements assuming all of the particles as being outgoing. Energy-momentum 
conservation will then force the energy of some of the particles to be negative. 

Notice that the following equations hold for generic chiral spinors (not necessarily solutions of 
a Dirac equation) which satisfy Ea. lTPJ) : fOjl . (1A~TU1) . (lA~T3i (lA~T4l) . (lA~T5l) . (lA~T6l) . 

The polarizations for vectors with momentum p, as defined in the text: 



e ± (p, k) -7 

enjoy the following properties: 

ej(p,*0 = ( 



>±|7„|fc±) 
V2(kT\p±) 



V2 



IpTX&Tl + |*±>(p±|), 



e (p, 



e (p, fc 



e (p, k 



e ± (p, k 



p = e ± (p,k) ■ k = 0, 



■e ± (p,k') = 0, 



e*{p,k') = -1, 



e ± (p',fc) = 0, 



<rF(Jfe,jfe') = 0, 



P/Lt^y + Pl/kf! 

p ■ k 



(A.17) 
(A.18) 

(A.19) 
(A.20) 
(A.21) 
(A.22) 
(A.23) 
(A.24) 

(A.25) 
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B Appendix: A QED Example 



In this Appendix we will describe as a simple application of the Helicity Amplitude technique the 
process of electron-positron annihilation into a photon pair. Two diagrams contribute to the process 
- t-channel and u-channel fermion exchange (see Figure 12). If q, q are the momenta of the incoming 
electron and positron, p 12 are the momenta of the two outgoing photons and h\ j2 are their helicities, 
the contributions of the two diagrams are the following: 

if? 

M * = T ^(q±\e h ^p 1 ,k 1 )(q-p 1 )e h ^p 2 ,k 2 )\q±), (B.l) 

[q-piY 

M « = T ^(q±\e h ^p 2 ,k 2 )(l-p 2 )e h ^p 1 ,k 1 )\q±). (B.2) 

{Q-P2) 2 

It is straightforward to check, using Eq. (jA.18|) . that if both photons have the same helicity then 
M = M u + Mt = 0. In fact by choosing the reference momenta to be equal to q (q) when the 
common photon helicities are the same as (opposite to) the electron helicity, we easily find that 
both M t and M t identically vanish. Therefore the only processes which contribute have photons 
with opposite helicity. By using the above criterion for the assignment of the reference momentum, 
and by assuming for definiteness that photon 1 has the same helicity as the electron, we can easily 
transform Eqs. (jB.l|) and (jB.2|) into the following expressions: 



-ze 



2 



[±)y/2 (=f)V2 



Mt = W V ,w-T =\ (v±\PiT)(qT\(<j-Px)\qT)(p2T\q±) (B.3) 



2ic 2 (pi ±\QT)(qT\p2±) = 2e 2 e ^ (m) 
(piT\q±)(q±\p2T) \(qp2) 

M u = 7 { ie \ {± 2f\\ l {T }f\\ (g±|g*>fa T| Cq-P2) \piT)(qT\q±) = 0. (B.4) 
2(gp 2 ) {qT\Pi±} (g±|pzT> 

Notice that even though in this gauge only the t-channel Feynman diagram is different from zero, 
it nevertheless contributes to the w-pole because of the singularity present in the definition of the 
polarization vector. This is a common feature of gauges defined in terms of arbitrary vectors, such 
as axial gauges. 

Squaring and summing over the various non vanishing helicity configurations, we finally obtain: 

£|M(e+e-- 77 )| 2 = 8e 4 ^±^. (B.5) 

pol Ul 
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C Appendix: Feynman Rules and the SU(N) Algebra 



Our Feynman rules follow the Bjorken and Drell conventions. In this Appendix we will collect them 
in the form which is most appropriate for their use with the helicity amplitude and the dual color 
expansion. 

By convention we will always assume all of the particles as outgoing and we will order the indices 
clock-wise. Positive- and negative-helicity quarks are represented by bra spinors: 

(?±|, (C.l) 

while positive- and negative antiquarks are represented by ket spinors: 

\q T ). (C.2) 

Positive- and negative-helicity gluons are given by the the positive- and negative polarization vectors 
introduced in the text and in Appendix A. 

The fermion and vector propagators are given, respectively, by: 

, -- 2 g, v 5 a \ (C.3) 
q 2 p 2 

where the indices i,j and a, b are the color indices of the fermion and adjoint representations, 
respectively. Since we will always be calculating helicity amplitudes, therefore dealing with physical 
external gluons, the choice of Feynman gauge is technically equivalent to any other choice. 

The fermion-antifermion-gluon vertex is given by: 

^ 9 (A°)y 7, (C4) 

The a/2 is a consequence of our choice of normalization for the A matrices and their algebra: 

[A a , A 6 ] = if abc \ c , tr(A a A b ) = 5 ab (C.5) 

(^ a )ij(^ a )ki = Sii5 jk - —SijSki- (C.6) 

a 

This is not the usual normalization, but this choice prevents the proliferation of powers of \/2 
which would otherwise appear in the calculation of the matrix elements and, independently, in the 
squaring of the color structures. All of these factors of 2 eventually cancel out, and our convention 
enforces these cancellations since the start. The three-gluon vertex is given by: 

^(Pl,P2,P3W 2M3 = 9 F(P1,P2,P3)km 3 , ( C - 7 ) 

+ (P2 - P3) ]. (C.8) 
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The \f2 is a consequence of our normalization of the SU (N) algebra. The three-gluon vertex breaks 
up into two pieces, corresponding to the two possible color structures this vertex can contribute to: 

g [tr(A ai A a2 A« 3 ) - tr(A a3 A a2 A ai )] F( Pl ,p 2 ,p 3 ) m . (C.9) 

In the calculation of a dual amplitude, where the diagrams are determined by a specified ordering 
of the gluons, we should use only the color structure corresponding to the assigned ordering: 

g tr(A ai A a2 A as ) F( Pl ,p 2 ,p 3 )^ 3 . (CIO) 

For example, the s-channel diagram contributing to the dual amplitude m(l, 2, 3, 4) would be given 
by: 

(73S) 2 E tr(X^X a -X b ) F(px,p 2 , P)p im v [~\ tr(\ a *\ a *\ b ) F(p 3 ,p 4 , P)^ v 
b r 



= {^ 2 gf tr(A ai A a2 A a3 A fl4 )[--^] F( Pl ,p 2 , P)^ 2 „ F(p 3 ,p 4 , P), 3 , 4 „, (C.ll) 

with P = pi + p 2 and where a term proportional to tr(A ai A a2 )tr(A a3 A a4 ) vanishes because of the 
anti-symmetry of the functions F. 

In an analogous way, and using the standard four-gluon vertex Feynman rule, one can write the 
four-gluon dual vertex, corresponding to the permutation (1,2,3,4): 

,|. tr (A^A a2 A a3 A a4 ) (2g^ 3 g^ 4 - g^ 4 g^ 3 - <W 2 <wJ (C.12) 

The sum of this vertex plus the other 5 obtained by permuting the indices gives the standard 
four-gluon vertex, as can be easily checked. 
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C.l Summary of Feynman Rules 



Here we summarize the Color truncated Feynman rules, where all the vertices are cyclically or- 
dered and all momenta are outgoing. Demonstrating that the subamplitudes m(gi, #2, #3, #4) and 
m(q, <7i, #2, q) are gauge invariant is an easy way to check the consistency of our conventions. 

• External, outgoing fermion, F, helicity ±: 

(F±|. (C.13) 

• External, outgoing anti-fermion, F, helicity ±: 

\FT). (C.14) 

• External, outgoing vector, momentum p, reference k, helicity ±: 

^ k ) = ^Mlt±r (ai5) 

^(P,k)-l = ± {k ^ p±) (\pT)(kT\ + \k±)(p±\). (C.16) 

• Fermion propagator, momentum q, in the direction of the fermion arrow: 

i f r (C.17) 



• Vector propagator, momentum p: 

— i 

p 

Fermion- vector-antifermion vertex, order (FVF): 



% ^f. (C.18) 



i-^ 7/1 . (C.19) 

Tri- Vector vertex, order (123), all momenta outgoing from vertex: 

i 71 tO* 1 ~ ^ 2 Www + {P2-Pz)^9^z + (P3-Pi) w ff w ]- (C.20) 

Quartic- Vector vertex, order (1234): 

9 2 
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D Appendix: Squares and Color Sums 



For the sake of definiteness we will choose the fermion color representation to be the fundamental 
representation. With our conventions the A matrices are hermitian and are normalized by 

tr(X a X b ) = 5 ab , (D.l) 

and satisfy the following identities which are useful in reducing the color sums of products of traces: 

jv 2 -i 1 

^ (A )i 2 j 2 = ^iij2^i'2ji ~ ~7r$iih^hh (D-2) 

0=1 

[tr(A ai A a2 ...A a ™)]* = tr(A a "...A a2 A ai ). (D.3) 

As a short-hand notation, we will define: 

( ai a 2 ...a n ) = tr(X ai X a2 ...X an ) (D.4) 

In squaring the gluon amplitudes and summing over colors, we will have to carry out sums of the 
following form: 

N 2 -l 

{a 1 a 2 ...a n ){b 1 b 2 ...b n )* (D.5) 

ai,...,On=l 

where {&} is a permutation of {a}. Using the cyclic property of the trace, the properties under 
complex conjugation and Equation (|C.6J) . we can always rearrange the permutation {b} so as to 
reduce the sum in the following way: 

N 2 -l N 2 -l 

nO"m n —\ ■ ■ ■ Q>m\) • Q"n—lQ"m n —i • ■ ■ Q"m\) 

ai,...,o n _i=l an oi,...,On_i=l 



- ^(oi...o t »_i)(fl mn _ 1 ...o mi )]. (D.6) 

The first term can have either of the following forms: 

^ (N 2 -l) 

(A^-iOn-iAa) = - — — — -(AiA 2 ), (D.7) 

a-n-l 

or: 

X; (Aia„-iA 2 a„-i) = (Ax)(A 2 ) - -^(AiA 2 ). (D.8) 

On-l 

In either case the sum over colors can be iterated using the formulas given so far, until all of the 
terms that will develop are proportional to {a\a{) = (N 2 — 1). The leading contributions in N come 
from {b} = {a}, in which case, to the leading order in N: 

N 2 -l 

X {aia 2 . . . a n )( ai a 2 . . . a n )* = N n ~ 2 (N 2 - 1) ( 1 + 0(1/N 2 ) ) (D.9) 

ai,...,a n =l 
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For the readers convenience we will collect here some useful formulas involving traces and squares 
of A matrices: 

(ab) = 5, 



ab 

N 2 - 1 

E(*°*% = — — 6y 
a 

J2(ab)(ab) = N 2 — 1 

ab 

( N 2 - 1 ~ 

Y,(abc){cba) = (N 2 — 2) 



abc 

y2(abc)(abc) 



'N 2 -l\ 



abc 



\ N 

J2f abc f abc = 2N(N 2 -1) 



abc 

As an explicit example we will now calculate the square of the 4-gluon amplitude: 

_/^f( 4 ) — V" (01020304)7^(1,2,3,4) = [01020304)777,(1,2,3,4), 

{2,3,4} {2,3,4}' 

where the prime indicates that only permutations inequivalent under reflection ( (1234) 
) should be considered, and where we introduced the following notation: 

[aia 2 ...a n ] = (aia 2 . . . a n ) + (-l) n (a n . . . a 2 ai), 

[aia 2 ...a n ]* = (-l) n [a„ . . . a 2 ai]. 
The sum over colors of the squared 4-gluon amplitude can be then written as: 

J2\M^\ 2 = J2 m(l,2,3,4)^[oi0 2 o 3 o 4 ] ([oi0 2 o 3 o 4 ]*m*(l,2,3,4) + 

col {2,3,4}' col 

[aia 3 a 2 a4]* m*(l, 3, 2, 4) + [0x020403]* m*(l, 2, 4, 3)) . 
It is simple to prove the following equations: 

^[01020304] [01020304]* = E[aia 2 a 3 a 4 ] 2 

col col 

^[01020304] [01020403]* = E[ aia 2 a 3 a4 ] 2 + iVfa^Os] 2 

col col 

^[01020304] [01030204]* = ^2[aia 2 a 3 a A ] 2 + A^[ai02a 3 ] 2 

col col 



Y,N[ ai a 2 a,} 2 = -2N 2 (N 2 -1). 

col 
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By using these equations and the Dual Ward Identity, Equation (|3.4j) . we therefore obtain: 
]T|M( 4 )| 2 = £ |m(l,2,3 l 4)| 2 (-X;^[oia2a8][aia a as]) 

col {2,3,4}' col 

= 2N 2 (N 2 -1) |m(l,2,3,4)| 2 

{2,3,4}' 

= iV 2 (iV 2 -l) J2 \ m ( l , 2 ^,A)\ 2 . (D-24) 

{2,3,4} 

Summing over the different helicity configurations, finally, gives Equation ()5.8j) . 

In analogous way one can show the vanishing of the subleading terms in the square of the 5-gluon 
process. As for the 6-gluon case, by using the DWI one can show that the square takes the following 
form: 

]T|M (6) | 2 = N\N 2 -l) (|™(123456)| 2 

col {23456} 

+ — m(123456)* x [m(135264) + m(153624) + m(136425)]) , (D.25) 

where the expressions of the dual amplitudes m(123456) for the various helicity configurations were 
given in Section 5 . 

Finally, we give here the structure of the amplitude squared for the processes {qqgg), {qqggg) 
and (qqgggg)- To keep the following formulae as simple as possible, we introduce the following 
notation for the quark sub-amplitudes: 

m(q, gi, gj,..., g L , q) = (I, J, . . . ,L), (D.26) 

where (J, J, ... , L) is an arbitrary permutation of (1, 2, . . . , 4). From the expansion of the amplitude 
in the usual color basis, 

M(q, gi ,...,g n ,q)= ]T (A : A 2 . . . X n ) l3 (1, 2, . . . , n), (D.27) 

{l,...,n} 

we obtain the following expression: 

( at 2 — 1) n ~ 1 

E \M(q, 9l ,...,g n ,q)\ 2 = £ ^ E ^(M, . . . ,n). (D.28) 

colors j=0 {l,...,n} 

For n = 2, 3, 4 the functions Hj are given by: 
• n = 2 

ff 1 (l,2) = |(l,2)| 2 (D.29) 
F (l,2) = -(1,2)*[(1,2) + (2,1)] (D.30) 
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• n — 4 



n = 3 

tf 2 (l,2,3) = |(1, 2, 3)| 2 (D.31) 

#i(l, 2, 3) = -(1, 2, 3)*[2 (1, 2, 3) + (1, 3, 2) + (2, 1, 3) - (3, 2, 1)] (D.32) 

#o(l,2,3) = (1,2,3)* E V,J,K) (D.33) 

# 3 (1,2,3,4) = 1(1,2,3, 4)| 2 , (D.34) 

H 2 (l, 2, 3, 4) = (1, 2, 3, 4)*[-3 (1, 2, 3, 4) - (1, 2, 4, 3) - (1, 3, 2, 4) 
-(2, 1, 3, 4) + (1, 4, 3, 2) + (3, 2, 1, 4) + (3, 4, 1, 2) 
+ (3, 4, 2, 1) + (4, 2, 3, 1) + (4, 3, 1,2)], (D.35) 

#x(l, 2, 3, 4) = (1, 2, 3, 4)* [M(l, 2, 3, 4) - M(4, 3, 2, 1)] (D.36) 
M(l, 2, 3, 4) = 3 (1, 2, 3, 4) + 2 (1, 2, 4, 3) + 2 (1, 3, 2, 4) + 2 (2, 1, 3, 4) 

+ (1, 3, 4, 2) + (1, 4, 2, 3) + (2, 1, 4, 3) + (2, 3, 1, 4) + (3, 1, 2, 4), (D.37) 

#o(l,2,3,4) = -(1,2,3,4)* ]T (I,J,K,L). (D.38) 

{I,J,K,L} 

The formulae for n = 2, 3 can be used to compare our results with the expressions already known. 
In doing this it is useful to apply the DWI to the functions Hj and use the gluino sub-amplitudes 
with non- adjacent fermions as auxiliary functions. 

The sub-amplitudes for the various helicity configurations were given in Section 5 . 



E Appendix: Numerical Evaluation of the Spinor Products 

To calculate the matrix element squared for a given process, it is frequently easier to evaluate the 
sub-amplitudes as complex numbers and then form the appropriate square using the color algebra 
of the previous Appendix. For all the processes discussed in this review, the sub-amplitudes can 
be calculated from sums of products of spinor products. Therefore, we need an algorithm for 
evaluating these spinor products. Given that both (ij) and [ij] are complex square roots of the 
Lorentz invariant SV,- = (pi + Pj) 2 , 

(u) = \f\Sij~\ exp(^), (E.l) 

M = ^\S~\ exp(ife)- (E.2) 

If both momenta having positive energy, the phase factor 0y is defined, in a popular representation 
of the gamma matrices, by 
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COS < 



I.) 



(pbj - 


pM) 






{pip~j - 


PjPt) 


JpTpj 


' <?■ • 



sin^- = ^ + + ^ ■ (E.3) 

Where p = (p° ± p 3 ) and since = 0, the spinor product for this representation of gamma 
matrices are undefined for a momentum vector in the minus 3 direction. If one or more of the 
momenta in (ij) have negative energy, 0y is calculated with minus the momenta with negative 
energy and then nir/2 is added to 0y where n is the number of negative momenta in the spinor 
product. The associated phase factor, 0^- , for [ij] can be calculated from Sy using the identity 
Sij = (ij) [ji]. The above identities can be used to evaluate the spinor products with approximately 
the same amount of computational effort as the evaluation of J . 

Similarly, the matrix element squared can usually be written simply as a sum of traces of 
momentum vectors, see If these traces are expanded, the resulting expressions are extremely 
cumbersome. However, the phase factors defined above can be used to evaluate these traces in an 
efficient manner. Consider the trace of a large string of light-like momentum vectors with all vectors 
having positive energy, then 

Tr(P 1 P 2 P 3 ---P 2n ) = [12](23)..-(2nl) + (12) [23] • • • [2nl] 



= 2 \j S 12 S 2 3 ■ ■ ■ S 2ni COS(012 - 032 + 034 - • • • - 01 2n)- (E.4) 

Where the identity for positive energy spinor products, 0jj = — 0ji, has been used. For traces 
involving 75, the corresponding identity is 

TriP 1 P 2 P 3 ---P 2nl5 ) = [12](23)---(2nl) - (12) [23] • • • [2nl] 



= -2i y S 12 S 23 . . . S 2nl Sin(0i2 - 032 + 034 - • • • - 01 2n)- (E.5) 

If m of the vectors in the string have negative energy, multiple these vectors by (—1) and use all 
the resulting positive energy vectors to calculate the above trace. The original trace is obtained by 
multiplying this answer by (— l) m . Traces involving vectors which are massive can be also treated 
by writing the massive vector as a sum of two light-like vectors (usually the decay products 
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Figure 1: The zero-slope limit of the four gluon string diagram in terms of Feynman diagrams. 




Figure 2: QED-type diagrams. 
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Figure 4: The color flows for quark-antiquark scattering. 
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Figure 5: The color flows for quark- ant iquark scattering with emission of two gluons. 
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Figure 8: The color ordered gluonic current. 
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Figure 9: A graphical representation for the Berends-Giele gluonic recursion relation. 
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